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Abstract. For a compact Riemannian manifold with boundary, endowed with 
a magnetic potential a, we consider the problem of restoring the metric g 
and the magnetic potential a from the values of the Marie action potential 
between boundary points and the associated linearized problem. We study 
simple magnetic systems. In this case, knowledge of the Mane action potential 
is equivalent to knowledge of the scattering relation on the boundary which 
maps a starting point and a direction of a magnetic geodesic into its end point 
and direction. This problem can only be solved up to an isometry and a gauge 
transformation of a. 

For the linearized problem, we show injectivity, up to the natural obstruc- 
tion, under explicit bounds on the curvature and on a. We also show injectivity 
and stability for g and a in a generic class Q including real analytic ones. 

For the nonlinear problem, we show rigidity for real analytic simple (g, a), 
rigidity for metrics in a given conformal class, and locally, near any (g, a) £ Q. 
We also show that simple magnetic systems on two-dimensional manifolds are 
always rigid. 
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1. Introduction 

1.1. Statement of the problem. Let M be a compact manifold with boundary, 
endowed with a Riemannian metric g, and let ir : TM — + M denote the canonical 
projection, tt : (x, £) i— » x, x £ M, £ £ T X M. 

Denote by ojq the canonical symplectic form on TM, which is obtained by pulling 
back the canonical symplectic form of T*M via the Riemannian metric. Let H : 
TM -> R be defined by 

H{v) = \\v\% v£TM. 

The Hamiltonian flow of H w.r.t. ojo gives rise to the geodesic flow of M. Let f2 be 
a closed 2-form on M and consider the new symplectic form ix> defined as 

The Hamiltonian flow of H w.r.t. u> gives rise to the magnetic (or twisted geodesic) 
flow ip* ■ TM — > TM. This flow models the motion of a unit charge of unit mass in 
a magnetic field whose Lorentz force Y : TM — > TM is the bundle map uniquely 
determined by 

n*(Z,v) = <y*(Z),v) (i.i) 

for all x £ M and £, 77 £ T X M. 
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Magnetic flows were first considered by V. I. Arnold in [4] and by D. V. Anosov 
and Y. G. Sinai in [3]. As shown in [3 EH EH E3 [HI [28] , they are closely related to 
other problems of classical mechanics, mathematical physics, symplectic geometry, 
and dynamical systems. 

It is not hard to show that the generator G M of the magnetic flow is 

where G is the generator of the geodesic flow, and that every trajectory of the 
magnetic flow is a curve of the form t ^ (7(t),7(i)), where 7 is a curve on M 
which satisfies the equation 

v^7 = r(7), (1.2) 

which is nothing but Newton's law of motion. Such a curve 7 is called a magnetic 
geodesic. Note that time is not reversible on the magnetic geodesies, unless £7 = 0. 
If £7 = we recover the ordinary geodesic flow and ordinary geodesies. 

When £7 is exact, i.e. Q = da for some magnetic potential a, the magnetic 
flow also arises as the Hamiltonian flow of H(x,p) = \{j> + a) 2 g with respect to the 
standard symplectic form of T*M. The function H is the symbol of the semiclassical 
magnetic Schrodinger operator. 

Since the magnetic flow preserves the level sets of the Hamiltonian function H, 
every magnetic geodesic has constant speed. 

Unlike the geodesic flow, where the flow is the same (up to time scale) on any 
energy levels, the magnetic flow depends essentially on the choice of the energy 
level. We fix the energy level H~ 1 (l/2), thus considering the magnetic flow on 
the unit sphere bundle SM of M , in consequence we consider only the unit speed 
magnetic geodesies. Note that fixing the energy level to be SM is no restriction at 
all, since one can obtain the behavior in any energy level by considering the flow 
on SM upon changing £7 by AO, where A 6 R. 

We define the action A(x, y) between boundary points as a minimizer of the 
appropriate action functional, see (jl.6p and Appendix A. In the case £7 = 0, A(x, y) 
coincides with the boundary distance function d g (x,y). In this case, we cannot 
recover g from d g up to isometry, unless some additional assumptions are imposed 
on g, see, e.g., [9]. One such assumption is the simplicity of the metric, see, e.g., 
[22" 1 l33 l l37 l [38], We consider below the analog of simplicity for magnetic systems. 

Let A stand for the second fundamental form of dM and v{x) for the inward 
unit normal to dM at x. We say that dM is strictly magnetic convex if 

A(x,Q>(Y x {t),v(x)) (1-3) 

for all (x, £) G S(dM) (see Appendix A for an explanation). Note that if we replace 
6 by — £, we can put an absolute value in the right-hand side of (| 1 . 3[) . In particular, 
magnetic convexity is stronger than Riemannian one. 

For x G M, we define the magnetic exponential map at x to be the partial map 
exp£ : T X M — > M given by 

exp£(t£) = t>0,£eS x M. 

It is not hard to show that, for every x € M, exp^f is a C 1 -smooth partial map on 
T X M which is C°°-smooth on T X M \ {0} (see Appendix A). 
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Definition 1.1. We say that M is simple (w.r.t. (<?, f2)) if dM is strictly magnetic 
convex and the magnetic exponential map exp£ : (exp^) _1 (M) — *• M is a diffco- 
morphism for every x 6 M (cf. the definition of a simple Riemannian manifold 
[22]). 

In this case, M is diffeomorphic to the unit ball of W 1 (so we can assume that 
M is this ball); therefore, £1 is exact, and we let a be a magnetic potential, i.e., a 
is a 1-form on M such that 

da = fl. (1.4) 

Henceforth we call (g,a) a simple magnetic system on M. We will also say that 
(M,g, a) is a simple magnetic system. 
Given x, y 6 M, let 

C(x, y) = {7 : [0, T] -> M : T > 0, 7(0) = x, 7(T) = y, 7 is absolutely continuous}. 
The time-free action of a curve 7 £ C(x, y) w.r.t. (g, a) is defined as 

A( 7 ) = A 3 , Q ( 7 ) = i^ | 7 ( t )|2 dt+ I T _y a . (i.5) 

For a simple magnetic system, unit speed magnetic geodesies minimize the time- 
free action (Lemma IA. 51 in Appendix A). More precisely, 



A(x,y):= inf A( 7 ) = A( 7X)2/ ) = T x , v - 

■y£C(x,y) 




where ^ x . y : [0, T XtV \ — > M is the unique unit speed magnetic geodesic from a; to y. 

The function A(x,y) is referred to as Mane's action potential (of energy 1/2), 
and we call the restriction A|a^/ x aM the boundary action function. 

Now, we state the boundary rigidity problem in the presence of a magnetic field 
as follows: To which extent is a magnetic system (g, a) determined by the boundary 
action function? 

To be more precise, we say that two simple magnetic systems {g, a) and (<?', a') 
are gauge equivalent if there are a diffeomorphism / : M — > ill, which is the identity 
on the boundary, and a function tp : M — > R, vanishing on the boundary, such that 
.9' = f*9 an d oi 1 — f*a + dip. Observe that gauge equivalent magnetic systems have 
the same boundary action function. 

Now, we rephrase the above problem as follows: Given a simple magnetic system, 
is any other simple magnetic system on the same manifold gauge equivalent to the 
former as soon as it has the same boundary action function? If so, we call the given 
magnetic system magnetic boundary rigid. 

Surely, this problem can be considered under various natural restrictions. For 
example, we can consider it for a fixed Riemannian metric and try to restore a mag- 
netic potential, or, vice versa, fix a magnetic potential and try to restore a metric, 
or consider the problem for metrics in a fixed conformal class, etc. In particular, 
for the zero magnetic potentials we recover the ordinary boundary rigidity problem 
for Riemannian metrics (see recent surveys of the latter in [101 IMj ) ■ 

For simple magnetic systems, knowledge of the action A(x, y) between boundary 
points is equivalent to knowing the scattering relation, see section 12.21 For non- 
simple systems, the problem of recovering (g, f2) from the scattering relation is the 
natural problem to study. The scattering relation appears naturally in the study 
of the scattering operator in K" with g Euclidean outside a ball, and a compactly 
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supported. Namely, for non-trapping metrics, the scattering operator associated 
to the semi-classical magnetic Schrodinger operator is a Fourier integral operator 
with canonical relation that determines the scattering relation on a large sphere, 
[H [15] . (It should be noted that for magnetic Schrodinger operators the resolvent 
is also a Fourier integral operator.) 



1.2. Description of the results. In section [2j we show that for simple magnetic 
systems, the action determines the jets of g and a in boundary normal coordinates. 
We define the scattering relation and show that for simple magnetic systems, it 
determines A(x, y) on the boundary, and vice-versa. We also show that one can 
recover the volume from A.{x, y). 

In section 02 we study the linearized problem. This reduces to the magnetic ray 
transform /. We show that potential pairs (see Definition 13. 5p belong to the kernel 
of /. We say that / is s-injective, if the kernel of / consists only of potential pairs. 

In section 21 we show that the normal operator N = 1*1 is a pseudo-differential 
operator in the interior of M, elliptic on solenoidal pairs that are an orthogonal com- 
plement of the potential pairs. We construct a parametrix of ./V; near the boundary, 
additional arguments are needed. This parametrix recovers the solenoidal projec- 
tion f s given Nf, up to a smoothing term. We show that s-injectivity implies a 
stability estimate, uniform near any (g,a), in appropriate spaces, see Theorem l4.3l 
We consider in this section and in Appendix B real analytic (<?, a) and show that 
then / is s-injective. A crucial element of the proof is that N is an analytic pseudo- 
differential operator in the interior of M. This is delicate since the magnetic ex- 
ponential map is only C 1 smooth when £1 ^ 0, even in the analytic case. To 
handle this, the analysis is done in polar coordinates. The s-injectivity for real 
analytic magnetic systems and the uniform stability estimate imply s-injectivity of 
the magnetic ray transform for generic (g,a). 

In section we show that / is s-injective for simple magnetic systems with an 
explicit bound on the curvature and £1. This relies on an analog of Pestov's identity 
for our case, see [30l [33] that goes back to [231 El], see also the references there. 

In section [SI we consider the non-linear magnetic rigidity problem. We prove 
rigidity in a given conformal class and rigidity within real analytic systems. We 
also show that if a simple Riemannian manifold (M, g) is boundary rigid (within 
the class of Riemannian metrics), then it is also magnetic boundary rigid. In this 
section we also study the local problem. We show rigidity near any (g,a) in the 
generic class Q using the analysis of the linear problem. This does not directly follow 
from the implicit function theorem, and the stability estimate in Theorem 14 . 31 plavs 
a crucial role. There is an essential difficulty compared to the Riemannian case 
a = 0, [37, 38J, since we cannot decouple g and a in the linearization argument. 
This difficulty is resolved by Lemma [6771 

Section [7] is devoted to two-dimensional systems. Here we prove that two- 
dimensional simple magnetic systems are magnetic boundary rigidity. This gen- 
eralizes the boundary rigidity theorem of [32] . Our proof essentially resembles that 
in 32J, establishing a connection between the scattering relation of a magnetic sys- 
tem and the Dirichlet-to-Neumann map of the Laplace-Beltrami operator of the 
underlying Riemannian manifold. 
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2. Boundary determination, scattering relation, and volume 

2.1. Boundary jets of the metric and magnetic potential. Here we show 
that up to gauge equivalence the boundary action function completely determines 
the Riemannian metric and magnetic potential on the boundary of the manifold 
under study. 

Lemma 2.1. If (g,a) and (g' , a') are simple magnetic systems on M with the same 
boundary action function, then 



where i : dM — * M is the embedding map. 

Proof. Given x € dM and £ € T x (dM), let t(s), — e < s < e, be a curve on dM 
with t(0) — x and f (0) = £. It is easy to see that 



Now, we prove that the boundary action function determines the full jets of the 
metric and magnetic potential on the boundary. This generalizes the corresponding 
results of [HHD]. 

Theorem 2.2. If {g,a) and (g',a') are simple magnetic systems on M with the 
same boundary action function, then {g' , a') is gauge equivalent to some (g, a) such 
that in any local coordinate system we have d m g\gM = d m g\gM and d m a\gM = 
d m a\oM for every multi-index m. 

Proof. Denote by v the inward unit normal to dM w.r.t. g. The "usual" boundary 
exponential map exp aM (p, t) — exjp p (tv(p)), p £ dM, t > 0, takes a sufficiently 
small neighborhood of the set dM x {0} in dM x M + diffcomorphically onto some 
neighborhood of dM in M. 

Let v' and exp' denote the corresponding objects for the metric g' . 

The map exp aM o (expg M ) _1 is well defined in some neighborhood of dM in M 
and is the identity when restricted to dM. We extend this map from a neighborhood 
of dM in ¥ to a diffeomorphism / : M — > M and put g — f*g' and a" — 
f*a' . By Lemma |2. 11 g and g induce the same Riemannian metric on dM and, by 
construction, at each point of dM the inward unit normal w.r.t. g coincides with 
the one w.r.t. g. Therefore, g\dM = g\dM- 

Next, by Lemma [2TTl a and a" induce the same 1-form on dM. Applying Lemma 
2.2 of [3S] to the form u> — a — a" we find a function ip e C^°(M) such that ui — dip 
induces the zero form on every sufficiently short geodesic of g starting from dM in 
the normal direction. We now put a — a" + dip. 

Let us prove the equality of derivatives on the boundary. We will use the same 
argument as in [20 . We fix xq £ dM and introduce boundary normal coordinates 
(x',x n ) w.r.t. g near xq. By construction, the same coordinates are boundary 



i" g = i* g' , i*a — i*a , 



(2.1) 



inn — - — = ia-«(o- 




□ 
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normal coordinates w.r.t. g. Thus, the line elements ds 2 and ds 2 of these metrics 
are given by 

ds 2 = g LK dx[dx' K + dx^, 

ds 2 = g bK dx[dx' K + dx 2 n , 

where i, k vary from 1 to n — 1. Therefore, for h — g — g we have hi n = for 
i = 1, . . . , n. Also, by the construction of a, for f3 = a — a we have f3 n = 0. 
It now suffices to prove that 

d™h lK \ x = Xo =0, d™fa\ x=xo =0 form = 0,l,...; l,k= l,...,n-l. (2.2) 

The case of to = is granted. Assume that there is a least to > 1 such that 
(|2.2[) is not true. The Taylor expansion of h and (3 then implies that there is a unit 
vector £o G T Xo (dM) such that 

1 



-h LK (x)et; K - fr(x)e >0 (or <0) 



(2.3) 



for x n > and x' both sufficiently small and £ close to £o- (Here t and k vary from 
1 to n — 1 because = h n \ = and f3 n = 0.) 

For arbitrary x,y £ dM, let 7 = y x , y : [0, T] -> M (7 = 7^ : [0, f ] -> M) be 
the unit speed magnetic geodesic of the system (<?, a) (system (g, a)) from x to y. 

On the one hand, since 7 minimizes the time-free action w.r.t. (<?, a), we have 



9Mt))Y(t)f(t)dt + ±T- 



a i (rf(t)y/ i (t)dt. 



On the other hand, since 7 minimizes the time-free action w.r.t. (g,a), 



Hx,y) = - 



Therefore, 







[9v(l(t)W(t)j\t)dt + -T- 



i%( 7 (t))7 l Wy'W-A(7(t))7 l W 



Similarly, we derive the inequality 
rf r 







^ y -(7(t))fW7 J W-ft(7W)fW 



dt < 0. 



dt > 0. 



(2.4) 



(2.5) 



Continuing the proof of the theorem, we now choose x = Xq and choose y = 5(s), 
where 6 : (— e,e) — > dM is a smooth curve with 6(0) — x$ and 6(0) — £o- Then 
for s > sufficiently small we see that, in view of (|2.3p . we cannot simultaneously 
have l|2.4p and (|2.5[) . This contradiction concludes the proof of the theorem. □ 



2.2. Scattering relation. Now, we define a scattering relation and restate our 
problem in terms of the scattering relation. 

For (£,£) G SM, let j x ^ : [£~ (x , £) , l(x , £)] -»Mbea magnetic geodesic such 
that 7^(0) = x, 7^(0)'= £, and *y x£ (£-{x, £)), 7sjC (*(ar, £)) £ 9M. Clearly 
the functions £~(x,£) and €(cc,£) are continuous and, on using the implicit function 
theorem, they are easily seen to be smooth near a point (x, £) such that the magnetic 
geodesic ^^(t) meets dM transversally at t — £~(x,£) and t = £(x,t;) respectively. 
By (|1.3p and Lemma IA.6I in Appendix A, the last condition holds everywhere on 
SM \ S(dM). Thus, £~ and I are smooth on SM \ S(dM). 
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Let d+SM and d-SM denote the bundles of inward and outward unit vectors 
over dM: 

d+SM = {(x,£) G SM : x G dM, (£,v(x)) > 0}, 
d-SM = {(x, £) G SM : x G dM, (£, i/(x)> < 0}, 

where v is the inward unit normal to dM. Note that d(SM) = d+SM U d-SM 
and d + SM n 0_0Jlf = S(<9M). 

Lemma 2.3 (cf. 34, Lemmas 3.2.1, 3.2.2]). For a simple magnetic system, the 
function L : d(SM) — > R, defined by 



t(x,0 if (x,t)ed+SM, 
r(x,0 if (x,0edSM, 



is smooth. In particular, £ : d+SM — > R is smooth. The ratio ^^fl) * s uniformly 
bounded on d(SM)\S(dM). 

Proof. Let p be a smooth nonnegative function on M such that dM = p _1 (0) 
and | grad p = 1 in some neighborhood of dM. Put h(x,£,t) = p{jx.^(t)) for 
G S(5M). Then 

/i(x,C,0) = 0, 

|£0r,£,o) = K»0,O> 

— (x, £, 0) = Hess, p(i, + Ms), ^(0>- 
Therefore, for some smooth function R(x,£,t), 

h(x, £, t) = (v{x), £)t + \ (Hcss x p{i, + (v(x),Y(t))) t 2 + R(x, £, i)i 3 . 
Since £, L(a;, £)) = 0, it follows that L = L(a;,£) is a solution of the equation 
F(x, £ L) := (1/(1), + i (Hess, £) + (i/(a:), £ + i2(a;, £, i)L 2 = 0. (2.6) 

By {OJ, for (x,0 £ S(dM) 

^-(x, £, 0) = \ (Hess, + nO)) 

= i(-A( a ; s O + Ma:),nO»<0. 

Now, the implicit function theorem yields smoothness of L(a;,£) in a neighborhood 
of S(dM). Since L is also smooth on d(SM)\S(dM), we conclude that L is smooth 
on d(SM). 

Next, from JMD we get for (a;, £) G <9(SM) \ S(9M) 



l -(Kess x + (i/(a;), + iZfo £, t)L(i, 



Again by (| 1 . 3[) this yields boundedness of the ratio i^ x f\^ f° r ( x i0 sufficiently 
close to S(dM) (where L is sufficiently small), and clearly implies boundedness of 
the ratio on the whole set d{SM) \ S{dM). □ 
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Definition 2.4. The scattering relation S : d+SM — > d-SM of a magnetic system 
(M, g, a) is defined as follows: 

S(x,0 = (7x,f(*(s,0) s 7x,c(^0))- 

The restricted scattering relation s : d+SM ^ M is defined to be the post- 
composition of the scattering relation with the natural projection of d-SM to M, 
i.e., 

s(x,£)=7* )S Mz,0)- 
By the preceding lemma, 5 and s are smooth maps. 
The next lemma generalizes the well-known assertion of [22] . 

Lemma 2.5. Suppose that (<?, a) and (</, a') are simple magnetic systems on M 
such that g\dM = g'\dM- If the boundary action functions A\qmx6m and A'\qmxOM 
of both the systems coincide, then the scattering relations S and S' of these systems 
coincide, S = S' . 

In the opposite direction, we have the following: 

Lemma 2.6. Suppose that (g,a) and (g' ,a') are simple magnetic systems on M 
such that g\dM = g'\dM and i*a = i*a! . If the restricted scattering relations 5 
and s' of the systems coincide, s = s' , then the boundary action functions of these 
systems coincide, A\ dM xdM = A'\ dM xdM- 

To prove these lemmas, we need one fruitful result. 

Lemma 2.7. If {M,g,a) is a simple magnetic system, then for x,y G dM 
dA(x,y) 



dA(x,y) 



= -(7*,i/(0),e> + «(0 for £ G T x {dM), 
= (%, y (T),ri)-a(ri) for 77 G T y (dM), 



dn 

where j x , y : [0, T] — ► M is the unit speed magnetic geodesic from x to y. 

Proof. Fix x, y G dM, x ^ y, and £ G T x (dM). Let r(s), — e < s < e, be a curve on 
(9M with t(0) = 2; and f (0) = £. For every s, put 7(i, s) = 7 T (s),y(i) 5 < t < T s , 
denoting 7 = j(t, 0) and T = To, and consider 

c(t,s)= 1 {^t,s), 0<t<T. 

Each curve c(-,s) is defined on the interval [0,T] and its length is exactly T s . We 
have 

dA(x,y)_dT Ho) d( f 



7a 



3£ ds ds 

Using the first variation formula for length and the fact that c(t,0) = 7(f), we 



have 

AT 

L (0) = -<7(0),0- / (V 7 7,W 



or. 

ds 



where V(t) — ff (i, 0) is the variation field of c(t, s). Using the equation of magnetic 
geodesies (|1.2p and the definition (jl.ip of F, we obtain 

-r(o) = -(7(o),0- / n( 7> 7) 
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On the other hand, it is easy to see that 



n}= H- / Q(V,i)dt. 



d 
ds 

This gives the first formula of the lemma. A similar calculation gives the second 
formula. □ 

Proof of Lemma \2.5\ By Lemma |2.1[ a(v) — a'(v) for all v G T(dM). Using 
Lemma [2771 and the fact that both metrics are the same on the boundary, we easily 
conclude that the scattering relations are the same. □ 

Proof of LemmaWM, Take x G dM and define s x : d+SM \ S{dM) -> dM \ {x} 
by &x{Q = This map is a diffeomorphism. Consider its inverse s" 1 : 

dM \ {x} — * d + SM. By Lemma l2~7l for every y G dM, y ^ x, we have for all 
£ G T a (0M) 

^|^ = -fe 1 (y),0 9 +«(0. (2-7) 

The assumptions of the lemma imply that the right-hand side of (|2.7[) is the same 
for the second magnetic system; therefore, so does the left-hand side. Now, the 
claim of the lemma is immediate. □ 



Thus, for simple magnetic systems, the boundary rigidity problem is essentially 
equivalent to the problem of restoring a Riemannian metric and a magnetic poten- 
tial from the restricted scattering relation. 



2.3. Determination of volume. Here we show that the boundary action function 
determines the volume of the manifold. This generalizes the well-known assertion 
of [H Proposition 2.13]. 

Theorem 2.8. // (g, a) and (g' , a 1 ) are simple magnetic systems on M with the 
same boundary action function, then the volume Yo\ g M of M w.r.t. g equals the 
volume Voh,' M of M w.r.t. g' . 

Proof. By Santalo's formula (see (|A.4[) in Appendix A) we have 
Vo\ g M=—f l{x,Z)dn(x,£). 

Wn-l J d+SM 

Using 

/ a{x,S)dY? n -\x^)=Q (2.8) 

J SM 

and Santalo's formula again, we obtain 

Vol 9 M = [ A( 7a!iC ) d(i(x, £). (2.9) 

Wn-l Jd+SM 

In view of Theorem 12 . 21 and Lemma |2.5( we may assume that the right-hand side 
of this equality is the same for the magnetic system (</, a'). This yields the sought 
equality of volumes. □ 
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3. Magnetic ray transform 



3.1. Derivatives of the action function. Let (g, a) be a simple magnetic system 
on M. It is easy to see that there is an e > so that every magnetic system 
(g + h,a + (3), satisfying 

\\h\\c*<e, W\W<e, (3.1) 

is simple. 

Given h and (3 satisfying (|3.ip . consider the 1-parameter family (g s ,a s ) with 
g s — g + sh, a s =a + s/3, se [0,1]. 
Clearly, each of these systems is simple. 

Lemma 3.1. For x, y £ dM, 

dA g s ia s(x,y) _ 1 



ds 



where "f s is the unit speed magnetic geodesic from x to y w.r.t. (g s ,a s ). 
If 

h\dM = 0, /3\dM = 0, 

then 



d 2 A g s t0l »(x,y) 



ds 2 



<C(\\hf cl + \\/3f cl ), 



(3.2) 

(3.3) 
(3.4) 



with a constant C independent of h and (3 and C 2 locally uniform in (g,a). 
Proof. Define 



1 



(p(s,r) := A 9%QT (7 S ) = - / \j s (t)\ -r dt + -T s - 



1, 



Then 



dA g s , a s(x,y) dip dip 

dS = ^ (s ' s) + a7 (s ' s) 



(3.5) 



By Lemma IA.5[ unit speed magnetic geodesies minimize the time-free action; 
therefore, for a fixed r, A s t )Q ,t (<y s ) has a minimum at s = r, which yields 

dip 



ds 



(r,r) =0. 



(3.6) 



Next, 



#|7*W ) dt - I £« T = £ / M> ~ / P- ( 3 - 7 ) 



dr 



Combining ([3"3 ]) -([3"77 ]) gives (f3~2")l . 

Differentiating (|3.2p and using (|3.3p . we obtain 



d 2 A g s^s(x 1 y) 
ds 2 



1 



where ^ = d-f s /ds. Whence 
d 2 A g s^ a s(x,y) 



ds 2 



yY.htj)¥.'yi + hij¥.(^Y.¥.) 

- (V 7i A-)7i-A(V 7i 7i)} dt, 
<C(||/ l || cl + ||/3|| cl )(IKIIc + ||V 7 - 7s ||c). (3.8) 
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By the equation of magnetic geodesies, we have 

= Y(js), (3.9) 

s s 

where V stands for the covariant derivative related to g s , and Y is the Lorentz force 
associated with g s , a s . The initial conditions are given by 

ls (0)=x, 7.(0) = (exp^r 1 ^), (3.10) 

where exp^ stands for the magnetic exponential map associated with 17 s , a s . From 
(|3.9p and (|3.10|) we easily infer that 

IIV.IIC + ||V 7i 7.||c <C(||/i||ci + ||/9||cO- (3-11) 
Combining and (j3~TT|) leads to (|3TI|) . □ 

3.2. Magnetic ray transform. Let (M,g,a) be a simple magnetic system and 
<fi : SM — y R a smooth function on the unit sphere bundle. We define the magnetic 
ray transform of <fi to be the following function on the space of unit speed magnetic 
geodesies going from a boundary point to a boundary point: 

I# 7 )= [<f>:= [ T ^(t)^(t))dt, 



where 7 : [0, T] — y M is any unit speed magnetic geodesic such that 7(0) € dM and 
7(T) e <3M. Assuming that the magnetic geodesies are parametrized by d+SM , 
we obtain a map J : C°°(SM) -> C{d+SM), 

I<t>(x, = r (X ' ? <M^(x, 0) (a, 6 fl+SM. (3.12) 



In the space of real- valued functions on d+SM define the norm 

IH| 2 =/ tfdy. 

Jd+SM 

and the corresponding inner product. Here d[i(x,£) = (£, j/(.t)) cE 2 "~ 2 (see IA.4I in 
Appendix A). Denote the corresponding Hilbert space by L^(d + SM). 

Lemma 3.2. The operator I extends to a bounded operator 

I : L 2 (SM) -► Ll(d+SM). 

Proof. Indeed, it is easy to see that {1(f) 2 < CI(f> 2 , with some constant C indepen- 
dent of 6. Therefore, 



/ {I(t)) 2 d^<C I 1(f) 2 d^ = C I tfdY, 271 - 1 

Jd + SM Jd + SM JSM 

by the Santalo formula (|A.4j) . □ 
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3.3. Solenoidal and potential pairs. For a decomposition of a symmetric tensor 
field into a potential and a solenoidal part (relevant when Q — 0), we refer to [33"ll34j. 

In view of the linearization formula (|3.2p , we are mainly interested in / applied 
to functions 4> that are of the form 

feO^iWee+ftwe* (3.13) 

(and, more generally that are polynomials in £). Then, given a symmetric 2-tensor 
h and a 1-form /3, we set for (x, £) € d+SM 

l[h,0\(x,t)= f hiMt))nt)¥(t)dt+ [ T PMVWVdt, (3.14) 
Jo Jo 

where 7 = 7^, T = t{x,£). 

If F is a notation for a function space (C k , L p , H k , etc.), then we will denote 

by F(M) the corresponding space of pairs f = [h, 0], with h a symmetric covariant 

2-tensor and (3 a 1-form, and denote by !F(M) the corresponding space of pairs 

w = [v, (p\, with v a 1-form and (f a function on M. In particular, L 2 (M) is the 

space of square integrable pairs f = [h, /?], and we endow this space with the norm 

l|f|| 2 = J\\h\ 2 g + ^\^dVol, (3.15) 

with the corresponding inner product. (The choice of the factor (n — l)/2 will play 
its role in the proof of Theorem [5TH) In the space C 2 (M) we will consider the norm 

||w|| 2 = f {\v\ 2 g + v 2 ) dVol. (3.16) 

J M 

Clearly, the norm of a pair [h,0\ in L 2 (M) is equivalent to the norm of the 
corresponding polynomial (|3.13[) in L 2 (SM). Therefore, Lemma 15721 implies that 

/ : L 2 (M) ^ L*(d+SM) 

and it is bounded. 

Define Y :T*M -> T*M by 

Y(v) = -(Y?vj), i = 

Note that this definition agrees with the map Y : TM —* TM and the isomorphism 
between the tangent and cotangent bundles via the Riemannian metric. 

Clearly, I vanishes on functions <j>(x,£) = G^ipix^) if vanishes for x E dM. 
To find a class of such functions that are polynomials of £ of degree at most 2, 
assume that 

G^[ Vi (x)e + v (x)] = hijwee w. 

Since 

G,(v t (x)e + <p(x)) = eivijC + Vi Yj + Vti ] 

= (d s v) ij ee + (v, j -Y(v) j )? ) (3.17) 

we get 

h = d s v, P = dcp-Y(v), (3.18) 
where d s v is the symmetric differential of v. We used here the fact that the even 
part of (|3.17[) w.r.t. £ £ S X M determines the quadratic form (d s v)ij^ l ^> , while 
the odd part determines /3(£). Next, knowing (d s v )ij£ l £-' for n(n + l)/2 generic 
£ <E S X M is enough to recover v; similarly knowing /?(£) for n linear independent 
£ G S X M is enough to recover (3. We have the following stronger statement. 
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Lemma 3.3. Fix x G M . Given any open subset V of S X M , there exist N — 
n(n+l)/2+n vectors £ m € V, m = 1, . . . , N, such that [h, (3] is uniquely determined 
by the values oftp(x,£) := hij{x)^ l S, J + [3j{x)^ 3 at £ = £ m , m = 1, .. . , N. 

Proof. Since x is fixed, we denote V'(C) = x/j(x,£). We show first that ip(£,), known 
for all £ € S X M, determines [h, /?]. Since VKO is a linear functional of [h, 0], and the 
latter belongs to a linear space that can be identified with M. N , it is enough to show 
that = implies h = 0, j3 = 0. By replacing £ by — £, we get that hij^ 1 ^ = 0, 
/Sj^- 7 = 0. The second relation easily implies that /3 = 0. The first one implies h = 
easily as well (see, e.g., [34] , where also a sharp estimate is established). 

Next, assume that ip(£) = in V. Then ^(£) = on S X M by analytic contin- 
uation, thus [h, P] = 0. Finally, since for any fixed £, is a linear functional 
belonging to (TSL N )' = M. N , and {"0(C); £ £ ^} is a complete set, there exists a basis 
0(a), fc= l,...,iV, in it. □ 

Remark 3.4. We can say a bit more. Since the determination of [/i, /?] is done by 
inverting a linear transform, one can choose £ m continuously depending on a;, if x 
belongs to a small enough set X, such that the map {ij>(x, £, m (x)), m = 1, . . . , N} 
[h, [3] is invertible with a uniform bound on the inverse. Then this can be extended 
to compact sets X . 

Definition 3.5. We call a pair [h, 0] € L 2 (M) potential if the equations (|3.18| hold 
with [v,<p] e Hl(M). 

This can be written as follows: 

3)-*©; 

where d s stands for the symmetric derivative acting on covector fields and d is the 
usual differential acting on functions (which coincides with the symmetric differen- 
tial on functions). 

Clearly, potential pairs satisfy 

l(d[v,(p]) = 0- (3-19) 
This follows from (|3.17|) if v, ip are smooth, and follows by continuity for general v, 
tp, once we establish the mapping properties of N = 1*1 below. 

We will relate potential pairs to the non-linear problem. Let (g,a) and (g',a') 
be two gauge equivalent pairs, i.e., g' = f*g and a' = f*a + dip with some dif- 
feomorphism / : M — > M, fixing dM, and some function ip vanishing on dM. 
Linearize this near / = Id and ip = Q. In other words, let f T be a smooth fam- 
ily of such diffeomorphisms with fo = Id and let ip T be a smooth family of such 
functions with ip° — 0. Let g T — f*g, a T = f*a + dip T , and we will compute 
the derivatives at r = 0. It is well known [3~il (3.1.5)] and is easy to calculate 
that dg T /dT\ T= o = 2d s v, where v = df T /dr\ T= o. Let dip T /dr\ T= o = tp. Since 
a J = ( a j ° f^dfl/dx 1 + dip T /dx l , we get 



aj = v 3 djCti + ctjdiV 3 + ipj 

= v j a^ 3 + a.jV 3 ti + tp,i 

= v j a^ 3 + (a,V),i - v j a>j ti + 

= -{da)ijiP + (<x,V + ip) yi . 
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Using (|1.4|) . ([1.1)1 . and treating v as a 1-form (by lowering the index), we obtain 
This shows that 



d_ 





'1 r r 


T = 


2 5 '~ a . 



= d[v, — (a, v) — ip], where v = 

UT 



d^ 

' ^ = 1~ 
t=o dr 



T = 

(3.20) 

Since ip can be arbitrary (as long as it vanishes on dM) , we see that the linearization 
of {g 1 , a') near / = Id, ip = is given by d[v, ip] with v, ip that can be arbitrary. 
For future references, let us mention that by (|3. 17)) 

G M (w,e) = (dw,0 =► ^(w( 7 ), 7 ) = <dw( 7 ), 7 ) (3.21) 

for any unit speed magnetic geodesic 7 (t), where w = [u, ip], and we used the 
notation (f, £) = (h,£, 2 ) + where f = [ft,/3], and similarly (w, £) = («,£) + "y 9 - 

Definition 3.6. We say that I is s-injective if I[h,(3] = with [h,0] G L 2 (Af) 
implies that is potential, i.e., 

h = d s v, P = -Y(v)+dip (3.22) 

with G Hj(M). 

Let [h, 0] be orthogonal to all potential pairs. Then 

77 — 1 1 

(h,d s v) + —^—((3,-Y(v)+d<p) } dVol = 

1 M 

for all v,ip vanishing on dM. Then 

t Ush - ^Y(f3),v) + ^{S0)<p} d ™ = 0, 
where S is the divergence. Therefore, 

77 — 1 

Sh — Y{fi) = 0, 6(3 = 0. (3.23) 

Definition 3.7. We call a pair [h, (3} solenoidal if the equations (|3.23|) hold. 
This can be written as 

h\ „ P [5 -My' 



0, 5 : 



2 



v /3 y ' \o 5 

Then d = —S*. 

In terms of the operators 8, d, the solenoidal pairs are defined as the ones in 
Ker<5, and the orthogonal complement of Ker<5 is Rand, consisting of potential 
pairs. Next, we will describe the projections to solenoidal pairs. 

We will show first that the operator — 8d is an elliptic second order (and clearly, 
a formally self-adjoint) operator, acting on pairs w = [v, ip], where v is a 1-form 
and ip is a function on M. First, notice that 

(-<Sdw,w) = Itdwf, w e C °°(Af int ); (3.24) 

thus, — dd is a non-negative operator. Let cr p (P) stand for the principal symbol 
of P. Then (]3 . 24[) implies that for any such w and a fixed x 

(-<7 p (5d)w,w) = K(d)w)|| 2 , (3.25) 
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where the inner product and the norm are in the finite dimensional space to which w 
belongs, i.e., they are as in (|3. 16|) before the integration. We will show that — (T p (Sd) 
is in fact a positive matrix-valued symbol for £ ^ 0. Since it is homogeneous in 
the dual variable £, it is enough to show that — a p (8d)~w = implies w = for 
£ ^ 0. In view of (f3T25|) . if -a p (8d)w = 0, then cr p (d)w = 0. Then a p (d s )v = 
(see (j4.6jl ). £<p — 0, where [v,<p] = w. It is well known [MllM], and easy to see 
directly, that d s is elliptic, therefore v = 0. Next, = implies ip = as well. 

Note that the Dirichlet boundary conditions are coercive for Sd, because the 
latter is positive elliptic. We will show that the kernel and the cokernel of this 
elliptic problem are trivial. If w belongs to the kernel of Sd, and satisfies Dirichlet 
boundary conditions, it has to be smooth and then (I3.24| holds for it as well. 
Then dw = and w = on DM . This easily implies w = in M by integrating 
(|3.21| along magnetic geodesies connecting boundary and interior points, and using 
Lemma \'3. 31 Using standard arguments, one easily checks that the cokernel of Sd 
equipped with Dirichlet boundary conditions, is trivial as well. Indeed, fix u in the 
cokernel. Since u is orthogonal to Sdw for all w G Cg°(M mt ), we get that Sdu = 0. 
This in particular implies that u has a trace on <9M, and choosing w's with a dense 
set of normal derivatives on dM, we show that this trace vanishes. Therefore, u = 
by what we proved above. Denote by (<5d)u the Dirichlet realization of Sd on M. 
The arguments above show that (5d)o is an invertible self-adjoint operator. 

We define the following projections 

V = d(Sd) D 1 d, S = ld-P. (3.26) 

Then f s = Si is solenoidal, Vf = dw is potential, and we have the orthogonal 
decomposition 

f = f s + dw 

into the solenoidal and potential parts. 

Since / vanishes on 7 3 L 2 (M), the s-injectivity of I is then equivalent to the 
following: I is injective on <SL 2 (M). 

3.4. Adjoint of /. For a fixed simple (g, a), consider 

I:L 2 (M)^Ll(d + SM) 

and consider its dual 

I* :Ll(d + SM)^L 2 (M). 

We will now find an expression for I*. Let f = [h, (3] and 4>(x, £) be as in (|3 . 1 3[) . 
Let ip(x,£) G C(d + SM). Then 

(Jf,V)= / 1>(x,t)dn{x,Q 

Jd+SM 

x / hiA^ Jt))il Mil At) + 0A^ At)W„ At) dt. 



o 



hij (7»,e (*))7i,« (*)7^ (*) + A (T/z,t (*))7i, 4 («) 
By Santalo's formula (IA.4j) of Appendix A, we get 

(7f>)= / [hij(x)?? +0i(x)e]lt(x,Z)dE 2n - l (x,t), 



ISM 

where V>"(x,£) is defined as the function that is constant along the orbits of the 
magnetic flow and that equals ip(x,£) on d + SM. Let dcr x (£) be the measure on 
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S X M. Then 



ft (a) / £V"(z,0<MOdVolOr). 

M JS X M 

Therefore, see also [35] . 

I*i>=\ [ e^'ti.^.©,^ f £V(s,0<MOl- (3.27) 

3.5. Integral representation for the normal operator. Let M\ D M be an- 
other manifold with boundary (a domain in K 71 , actually) such that M] nt D M. 
Extend 5, a to Mi. Then (Mi, 5, a) is still simple if Mi is close enough to M. 
Choose xq € Mi \ M, and consider the map 

(exp^)- 1 :M 1 ^(exp^ o )- 1 (M 1 ). 

It is C°° away from x — xo, therefore it is a C°° diffeomorphism to its image, if we 
restrict it to any compact submanifold M1/2 of M{ nt . We can assume that Mi/ 2 
has a smooth boundary, and that M^f 2 D M. We denote M1/2 by Mi again, and 
the map above gives us global coordinates in a neighborhood of Mi that can be 
identified with a subset of M. n with a smooth boundary. So in this section, Mi is 
considered as a subset of R n . If g, a are analytic, and M is analytic, then we will 
also assume that 9Mi is analytic. 

Denote by I\ the magnetic ray transform on L 2 (Mi), 

h : L 2 (A/ X ) - LKd+SMx), 

and denote by I* its dual. Set 

N = Ilh. (3.28) 

Extending all tensors as on Mi, we consider L 2 (M) as a subspace of L 2 (Mi). 
In this way we may consider TV as the operator 

N : L 2 (M) -> L 2 (Mi). 

We then define s-injectivity of N as in Definition 13.61 

Lemma 3.8. N : L 2 (M) — > L 2 (Mi) is s-injective if and only if I is s-injective. 

Proof. Let N be s-injective. Assume that It = with f G L 2 (M). Then /if = 
on d+SMi as well, therefore Ni = 0, which implies that f is potential. 
Now, assume that I is s-injective and let Nf = 0. Then 

0=(iVf,f) L2(Ml) = ||/if|| 2 , (a+SMi) , 

therefore, If = as an element of L^(d+SM) as well, hence f is potential. □ 

We will find an integral representation of N and will show that it is a 'I'DO in 
a neighborhood of M, following [3711381136] . 

Using (|3.27p . and replacing £ by v, we get the following. 



Proposition 3.9. 



Nf = 



N 22 h + N n p, -{N l2 h + Nop) 

n — 1 



(3.29) 
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wh 



ere 



(7Vn/3) 1 (x) = / «' da x (v) / A(7^ (*))■)£,„(*) dt, 

JS x Mx J 

(N 21 /3f j ' (x) = [ w'V <M«) / ^(7^(*))^,,W 

JS x Mi J 



(3.30) 



(JW ' (a) = / u 1 ^ da x (v) / ^(7x,«(«))7i,t,(*)^,»(*) 

In each of the integrals above, split the t-integral into two parts: 7+ correspond- 
ing to t > 0, and corresponding to t < 0. In /+, make the change of variables 
V = Jx,v{t) = ex Px{t v )- Then tv = (exp^) _1 (y), and t,v are C°° functions of x, 
y away from the diagonal x = y. We treat /_ in a similar way by making the 
substitution t' = —t, v' = —v. To this end, note that the simplicity assumption 
implies that the map 

exp£>~(tv) = tto ir\-v), t > 0, v G S X M V (3.31) 

is also a C 1 diffcomorphism for every x G Mi. 

We analyze the Schwartz kernel of iV near the diagonal below, but at this point 
we just want to emphasize that it is smooth away from the diagonal and therefore 
N has the pseudolocal property. 



4. Analysis of the operator N. Generic s-injectivity 

In this section, we analyze N and prove s-injectivity for analytic (g,a), and for 
generic ones. 

Here we first assume that [g, a] G C fe (M). As a result, various objects related 
to g, a will have smoothness l(k) such that l{k) — > oo as k — > oo. To simplify 
the exposition, we will not try to estimate l(k) (actually, l{k) = k — ko with ko 
depending on n only). We will say that a given function (or tensor field) / is 
smooth, or that / G C l , if such an I = l(k) exists, and I may vary from line to line. 
In particular, if k = oo, then I = oo. 

At the end of the section we analyze the case in which (g, a) is an analytic 
magnetic system. 

4.1. More about the magnetic exponential map at the origin. In order 
to study the singularities of the kernel of N near x = y, we need more precise 
information about the exponential map at the origin. One of the interesting features 
of the magnetic problem is that the magnetic exponential map is not C 2 unless 
f2 = by Lemma I A. 71 On the other hand, in polar coordinates, it is a smooth map. 
We are therefore forced to work in polar coordinates. 

Consider y — exp^(to) = 7 X <;(£), where v G S x Mi, and t > are such that 
y G Mi. Recall that Mi now is a subdomain of R". We are interested in the 
behavior of y for small |i| near a fixed xq, therefore, we can assume that we work 
in U — {x; \y — xq\ < e} with < e « 1 such that U is strictly convex w.r.t. the 
magnetic geodesies as well as w.r.t. the Euclidean metric. 

By Lemma IA.7| the map tv = £ i— > y is not C 2 in general. On the other 
hand, the map (t, v) t— > y is smooth, and respectively analytic, if g is analytic, 
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and in fact extends as a smooth/analytic map for small negative t as well, by the 
formula y = ^ x ,v{t)- Similarly, the function m(t,v;x) — (~f x ,v{t) — x)/t has the 
same smoothness, therefore 

Jx,v(t) ~ x = tm(t, v; x), m(0,v;x)=v. (4-1) 

We introduce the new variables (r, ui) £ R X S X U by 

r = t\m(t, v; x)\ g , uj — m{t, v; x)/\m{t, v; x)\ g . (4-2) 

Then (r, ui) are polar coordinates for y — x = rui in which we allow r to be ne- 
gative. Clearly, (r, uj) are smooth/analytic at least for e small enough, if g is 
smooth/analytic. Consider the Jacobian of this change of variables 

, d(r, uj) 
J:=det-^. 

d(t,v) 

It is not hard to see that J\t=o — 1, therefore the map K. x S X U 3 (t, v) i— > (r, oj) 6 
K.x S X U is a local C l , respectively analytic, diffcomorphism from (a neighborhood 
of 0) x S X U to its image. We can decrease e if needed to ensure that it is a (global) 
diffeomorphism on its domain because then it is clearly injective. We denote the 
inverse functions by t — t(r,uj), v — v(t,lu). Note that in the (r, uj) variables 

t = r + 0(\r\), v = w + OQr\), %, v (f) = u + OQr\). (4.3) 

Another representation of the new coordinates can be given by writing 

Exp£(t,v) = J x , v (t). 

Then 

r = ,ig»( f ) |ExpS( t ,„) - .1. , » - «Xt> £*$:^ , 

and 

(t,v) = (Exp^ 1 ) -1 (z + ra) 
with the additional condition that r and t have the same sign (or are both zero). 



4.2. Principal symbol of N. Since we showed that the Schwartz kernel is smooth 
away from the diagonal, and we want to prove eventually that N is a ^DO, it is 
enough to study the restriction of N on a small enough set U as above. We analyze 
7V22, for example. Perform the change of variables (t, v) t—> (r,ui) in (|3.30|) , to get 

(A^22^) 2 3 (x) = I lv l (r, (jj; x)v-' (r, lu; x)hij (x + ruj)w z (r, uj; x)w 3 (r, uj; x) 

x J _1 (r,w;a;) dr da x (uj), (4.4) 

where w(r, ui, x) = % tV (t), and supp/i C U. 

This type of integral operators is studied in Appendix [B] Applying Lemma IB. II 
and the remark after it to A^, see (|4.4I) . and proceeding similarly for the other 
operators N^i, we get 
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Proposition 4.1. Nki, k,l — 1,2, are ^DOs in M[ nt with principal symbols 



v P (N n y\x,Q 


= 2ir 




= 2tt 


<j p (N 21 y'^(x,o 


= 2vr 


a p (N 22 y'i'v(x,0 


= 2tt 



Sxll 



s*u 



(4.5) 



Here w • £ = t^Zj- The reason for a p (Ni 2 ) — 0, <r p (N 2 i) = is that they 
are integrals of odd functions. We therefore get that, as a of order — 1, 

cr p {N) — dia,g{a p {N 22 ), ^rj-o-p(iVii)). Note that the principal symbols of Nn, N 22 
are the same as in the case of "ordinary" geodesies, i.e., they do not depend on Y. 
Explicit formulas for a p (Nn), cr p (N 22 ) can be found in [371 [35] an d they are based 
on the analysis of the Euclidean case in [33] . Similarly, the principal symbols of d 
and S are independent of Y and are given by 

a p {d)[v,4>] = diag fy&j , a p (S)[h,0\ = (i'//,..- i,) . (4.6) 

It is well known J37J [38] , and follows immediately from Proposition 14. 1[ that N 22 
is elliptic on solenoidal tensors, i.e., a p (N 22 )[h, (3] = and a p (S)[h,0\ — imply 
[h,0] = 0. Similarly, Nu is elliptic on solenoidal (divergence free) 1-forms. As a 
result, N is elliptic on solenoidal pairs. 

We want to emphasize here that pairs solenoidal in M and extended as zero to 
Mi \ M (that we always assume), may fail to be solenoidal in M x due to possible 
jumps at dM . 

4.3. Parametrix of N. We proceed as in [37J[3S]. We will define the Hilbert space 
H 2 (Mi) as in [371 [38]. Let x — (x',x n ) be local coordinates in a neighborhood U 
of a point on dM such that x n = defines dM. Then we set 

„ n— 1 

n/iiW) = / (Ei^/i 2 + i^"/i 2 + i/i 2 ) dVo ^)- 

Ju i=i 

This can be extended to a small enough neighborhood V of dM contained in Mi . 
Then we set 

n 

ll/lli^Mi) = E W^fWHHv) + Wfllnw ( 4 - 7 ) 

3=1 

We also define the H 2 (Mi) space of symmetric 2-tensors and 1-forms, and also the 
H 2 (Mi) space of pairs f = [h,/3] which we denote by H 2 (Mi). Clearly, the latter 
is a Hilbert space and H 2 (A/i) C H 2 (Afi) C H^Mx). 

The space H 2 (Mi) has the property that for each f S H 1 (M) (extended as zero 
outside M), we have Nf £ H 2 (Aii). This is not true if we replace H 2 (Mi) by 
H 2 (Mx). 

Lemma 4.2. For any t — l,2,...,oo, there exists k > and a bounded linear 
operator 

Q : H 2 (Afi) — > SL 2 (M) (4.8) 
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such that 

QNf = t s + Kf for all f G H 1 (M), (4.9) 

where K : H 1 (M) -> «SH 1+t (M) extends to X : L 2 (M) -> SH'(M). If t = oo, 
then k = oo. Moreover, Q can be constructed so that K depends continuously on g 
in a small neighborhood of a fixed go 6 C k (M). 

Proof. We will first reconstruct the solenoidal projection fj^ from Nf modulo 
smooth terms. Recall that we extend f as zero outside M, and fjj- is the solenoidal 
projection of the so-extended f in M%. Next, we will reconstruct f s from f| f . We 
follow [371 [38]. 

As in [37[[38], we will work with ^PDOs with symbols of finite smoothness k 3> 1. 
All operations we are going to perform would require finitely many derivatives of the 
amplitude and finitely many seminorm estimates. In turn, this would be achieved 
if g G C k , Y G C k , k 3> 1, and the corresponding ^PDOs will depends continuously 
on g, Y . 

Since A" is elliptic on solenoidal pairs, we get that 

W := N + NoVm-l : C°°(M) -> C l (Mi) 

is an elliptic ^DO of order —1 in M mt , where Ao is any properly supported 
parametrix of (—A) 1 / 2 (having principal symbol 1 ). Next, Wi, restricted to 
a small neighborhood of dM\ in M\, is smooth because suppf C M. Therefore, 
there exists a left parametrix P to W such that PW - Id : L 2 (M) -> H*(Mi) 
with t ^> 1, if k ^> 1. Then Sm x {PW — Id)5M a has the same property, therefore 
Pi := Sm x P satisfies P\N = Sm x + K 2 , where K 2 has the smoothing properties 
above, therefore, 

P 1 m = f s Ml +K 2 i. (4.10) 

The next step is to compare i s and f^- . We have f s = + du, where u = 
wjkj — w, and the latter are the potentials related to f in Mi and M, respectively. 
Notice that d commutes with the extension as zero when applied to w because the 
latter vanishes on DM. Then u solves the boundary value problem 

(<5d)u = in M, u\ 9M = w Ml \ dM . (4.11) 

We need to express wm, \om in terms of ATf . Since f = outside M, relation (|4.10|) 
implies that 

- dw Ml = PyNf - K 2 f in M x \ M. (4.12) 

Let (2, £) G S(M[ nt \M int ) be such that the magnetic geodesic 7 x ,f (t), t > 0, can be 
extended in Mj nt \M int for t G [0, £1(2:, £)) and J x ,^{£i{x, £)) G dMr, and moreover, 
7 Xl { is transversal to dMi. Such (x,£) clearly exist if Mi is close enough to M 
and if 7^ is close to the outgoing magnetic geodesic normal to dM . Using (|3.21D . 
integrate ([4. 12|) along such j x ^ to get 

(w Ml (x),0= ({PiNi-K 2 i){ lx ^)^ x ^)dt. (4.13) 

Jo 

Here we denote (w, £) = Vj^ J + ip, where w = [v(x), <p(x)] with a 1-form v and 
a function <p, see also (|3.21l) . Similarly we define (f, £). By Lemma 13.31 for a 
fixed x, one can reconstruct v{x), tp from (w, £) known for finitely many £'s in any 
neighborhood of a fixed £, and this is done by inverting a matrix. Moreover, one 
can do this near any fixed x, and the norm of the solution operator is uniformly 
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bounded in x, see the remark after Lemma 13.31 By a compactness argument, one 
can construct an operator P 2 such that 

w Mi i x )\dM = p 2(PiN - K 2 )f. (4.14) 

Arguing as in [37], we see that 

P 2 Pi : H 2 (Afi) -> H 1/2 {dM) (4.15) 

is bounded. Let R : H t_1 / 2 (M) -> H*(M) be the solution operator u = Rh of the 
boundary value problem ddu = in M, u = h on dM. Then P~TU]) , (|4"TT|) . (l4~T4l . 
and (14. 15)) imply (see also [37]) 

f s = (Id + dRP 2 )PiNf + Kf, 

where K has the required smoothing properties. We apply S to the identity above 
and set Q := S(ld + dRP 2 )P 1 . 

To prove the last statement of the lemma, we note that all ^>DOs we work with 
depend continuously on g, a if k 3> 1. The same applies to 5, R and P 2 . □ 

4.4. Main results for C k coefficients. 

Theorem 4.3. Let (g, a) be a simple C magnetic system on M extended to a 
simple magnetic system on M\. Then, for k 3> 1, 

(a) Ker/n<SL 2 (M ) is finitely dimensional and included in C l (M), where I — > oo 
as k — > oo . 

(b) Assume that I is s-injective for the pair (g, a). Then 

with a constant C > that can be chosen to be uniform in (g, a) under a small 
enough C k perturbation. 

Proof. This theorem is an analog of [371 Theorem 2] except for the uniformity state- 
ment which is similar to an analogous result in [35]. Part (a) follows directly from 
Lemma |4~21 Part (b), without the last statement, can be deduced from Lemma |4~21 
as well, as in [37J . Finally, the proof of the statement about the uniformity of C is 
identical to that of [38j Theorem 2] . □ 

Without the assumption that I is s-injective, one has a hypoclliptic a priori 
estimate that can be obtained from (|4.16|) by adding C s ||/||h-=(a/i)7 > 0, to its 
right-hand side, see [37] Theorem 2(a)]. 

Corollary 4.4. The set of simple C k magnetic systems (<?, a) with s-injective mag- 
netic ray transform I is open in the C k topology, if k 3> 1 . 

The next lemma is a linear version of Theorem 12. 21 see also [38j Lemma 4] . 

Lemma 4.5. Let (g,a) be a simple C k magnetic system on M. Let f £ L 2 (M) be 
such that If = 0. Then there exists a C l pair w, with I — > oo as k — * oo, vanishing 
on dM , such that for f := f — dw we have 

d m f\ dM = 0, \m\ < I, (4.17) 

and if f = [h,f3], then in semigeodesic boundary normal coordinates, 

h in =p n = Q for all i. (4.18) 

If k = oo, then I = oo and, in particular, (|4.17|) holds for all m. 
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Proof. In view of Theorem 14. 3f a). we may assume without loss of generality that 

f e C l (M). 

We now show how to construct f satisfying (|4.18|) . We refer to [H] [35] for similar 
arguments. Let (x',x n ) be semigeodesic boundary normal coordinates. Then g ln = 
S in , P nn = Tf n = for all i. Let f = [h,P], w = [t>,<p]. Condition ([CPS)) is 
equivalent to 

Vi«„+V B «j =2htn, a„v-r> t = /3„. (4.19) 
The first system can be solved by setting first i = n and solving d n v n = /i nn , v n = 
for i n = by integration. Then we solve the remaining system of n — 1 ODEs of 
the form 

d n v L - 2T? n v K = 2h m ~ d,v n , i = 1, . . . ,n - 1, (4.20) 

with initial conditions v L = for x n = 0. Finally, we solve the second equation in 
(|4.19| with the same zero initial condition. This defines w near dM. To define w 
in the whole M, we replace w by xw, where x is an appropriate smooth cut-off 
function equal to 1 near dM and supported in a larger neighborhood of dM. 

To prove (|4.17p . we argue as in Theorem l2.2[ following [20]. Assume that (|4.17[) 
is not true at some xq G dM. Then by studying the Taylor expansion of (f, £) = 
hund* near x = xq, we see that there is £o G T(dM) such that (f, £) is either 
strictly negative or strictly positive for £ lying in some neighborhood of £o- This 
contradicts the fact that If = if we integrate over magnetic geodesies originating 
from xq with directions close to £o- D 

4.5. Analytic magnetic systems. Assume that M is an analytic manifold with 
smooth boundary dM that does not need to be analytic (that can always be 
achieved by choosing an analytic atlas; and in case we have a simple system, we can 
start with a fixed global coordinate system and do only analytic changes of vari- 
ables). We will show that then / is s-injective, if (g, a) are analytic. By "analytic", 
we mean real analytic, and we say that / is analytic in the set X, not necessarily 
open, if / is analytic in a neighborhood of X. Then we write / € A(X). Notice 
that one can construct an analytic Mi as before. 
The central result of this section is the following. 

Theorem 4.6. If(g 7 a) is an analytic magnetic system on M, then I is s-injective. 

We will give a proof at the end of this section. 

Theorem 4.7. Let g, a be analytic in M\. Then Nki, k,l ~ 1,2, are analytic 
^>DOs in Ml with principal symbols as in Proposition ^. 1\ 

Proof. Notice first that by (|3.30[) and the simplicity assumption, Nki have Schwartz 
kernels that are analytic away from the diagonal. Therefore, it is enough to prove 
the theorem for Ni^ restricted to an arbitrary small open subset of Mi. This, 
however, follows from (|4.4p and Lemma IB. 31 of Appendix 151 □ 

Lemma 4.8. Let g, a be analytic in M\ and assume that If = 0, f £ L 2 (M). 
Then f s G A(M). 

Proof. Consider the solenoidal projection f| fi = f dwMi of f (extended as 
outside M) on M x . Since 8f% h = and Nffa = in Mj nt , and since d and N 
together form an elliptic system of analytic "JDOs (we can apply an elliptic ^DO 
of order 2 to the left of N to make the new operator and 6 of the same order, 
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see also [371 13H]), we get that f| fi e A(M[ nt ). On the other hand, wjfj solves 
Mw^ = in M[ nt \ M, wjv/i =0 on dM\, and by elliptic boundary regularity, 
see |38[ Lemma 3] and references therein, we have that is analytic up to dM\, 
therefore f s Mi G A{Mx). 

Next, we have f s = if^ + du, where u = wm, — w, see also the proof of 
Lemma 14.21 Then u solves the boundary value problem (|4.1ip , and all we need is 
to prove that wmJsm is analytic. Note that in general, wji/j is not analytic across 
dM, because dw^ may have a jump there but it belongs to H 1 because d is 
elliptic; therefore the trace on dM is well defined. Since dwjvfj = fjvfj in M[ nt \ M, 
and it is analytic in its closure, in the notation of (|4.13ll . we have 

rti(x,£) 

Jo 

for x <E dM and £ in a small neighborhood of the unit exterior normal to dM. 
This, combined with Lemma 13.31 and the remark after it, shows that WM 1 \dM is 
analytic. □ 

Lemma 4.9. Under the assumptions of Lemma \4-5\ assume that g and a are 
analytic in a neighborhood of dM. Then f = in a neighborhood of dM. 

Proof. It is enough to notice that near dM, w in the proof of Lemma [4.51 is ob- 
tained by solving ODEs with analytic coefficients. Therefore, f is analytic in a 
neighborhood of dM, and it must vanish there by (|4.17|) . □ 

Proof of Theorem \4-b] We first find w as in Lemma 14.51 and Lemma 14.91 such that 
f = f s dw vanishes near dM, i.e., 

P = dw near dM. (4.21) 

Our next goal is to show that one can extend w = [v, <p] to the whole M analytically; 
then ()4.21j) would be preserved by analytic continuation (recall that f s is analytic 
by Lemma I4T51) . and this would imply f s = 0. 

Let u±(x, £) be the solution of the kinetic equation G^u — (f s , £) with u — on 
d±SM. Integrate to get 

Jo 

where £ + stands for I. Since H = 0, we have u+ +it_ = 0. For such that x is 

close enough to dM and £ close enough to — e n (in boundary normal coordinates), 
we have u+(x, £) = (dw, = (v(x),£)+(p(x) by (|4.21|) and (|3.21[) . This also implies 
— U-(x, £) = (v(x),£) + <p(x). If we replace ^ by — since — £ is close to e„, we get by 
(|4.2ip . — u~(x, — = —{v(x),£) + <p(x); therefore, u+(x, — ^) = —(v(x),£) + <p(x). 

In particular, for x close enough to dM and £ € S X M close enough to either e„ 
or -e„, 

<p(x) =u + , even (x,£), Vj(x) = d nj \r]\gU + ^ Ai {x,r]/\r]\g)\^ g=v (4.22) 

where u+, odd/even stands for the odd/even part of u + w.r.t. £. The derivative on 
the right-hand side of (14. 22[) can be written as Pu+. dd, where P is a first order 
differential operator on S X M with coefficients analytically depending on x and £. A 
direct differentiation shows that (Pu)j — d^u — ^j^dpu + £jU, and the first order 
part is clearly tangent to S X M. 
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We define y>(x, £) and v(x,£) by (|4.22|) on the whole SM. Since u+(x,£) is 
analytic on SM, so are ip(x,£) and u(a;, £). By (|4.22|) . in some open set they are 
independent of £, i.e., grad^(p(x, £) = 0, g£ad^v(x,^) = there, where grad^ stands 
for the gradient on S X M. Those equalities extend to the whole SM by analytic 
continuation; therefore, (|4.22|) defines ^-independent <p and v in the whole M, and 
they are analytic functions of x. □ 

We are ready now to state the generic s-injectivity result. 

Definition 4.10. For a fixed manifold M, we define Q k to be the set of simple C k 
pairs (g,a) with s-injective magnetic ray transform I = I g , a - 

Theorem 4.11. There exists ko > 0, such that for k > k$, the set Q k is open and 
dense in the set of all simple C k pairs (g, a) and contains all real analytic simple 
pairs. 

Proof. By Corollary 14. 4[ Q k is open. By Theorem 14.61 it is dense. □ 



5. Energy estimates method 

It is easy to see that if (M, g, a) is a simple smooth magnetic system and the 
magnetic ray transform of a smooth function <f> : SM — > R vanishes, then <f> is 
the flow derivative of a unique function u which is continuous on SM, smooth on 
SM \ S(dM), and vanishes on d{SM): 

G^u = cf), u\ 9(SM) = 0. (5.1) 

Indeed, u is defined by the formula 

u(x, = - / (x, 0) dt, (x,0e SM. (5.2) 

Jo 

In this section we will analyze the linear problem (|5.1[) for <f>(x,C) of degree at 
most 1 in £ (the linear problem for 1-tensors) and for <fi(x,£) of degree at most 2 in 
£ (the linear problem for 2-tensors). 

5.1. Semibasic tensor fields. We recall the notion of semibasic tensor field and 
its derivatives (we prefer to adhere to the notations of [12]). 

Let 7T : TM\{0} -> M be the natural projection, and let (3 r s M := n*r 7 s M denote 
the bundle of semibasic tensors of degree (r, s), where rJM is the bundle of tensors 
of degree (r, s) over M. Sections of the bundles (3 r s M are called semibasic tensor 
fields, and we denote the space of smooth sections by C°°(/?gM) (in particular, 
C°°(/3gM) = C°°{TM \ {0})). For such a field T, the coordinate representation 

T={T£;;£)(x,t) 

holds in the domain of a standard local coordinate system {x 1 ,^ 1 ) on TM \ {0} 
associated with a local coordinate system (x l ) in M. Under a change of a local 
coordinate system, the components of a semibasic tensor field are transformed by 
the same formula as those of an ordinary tensor field on M. 

Every "ordinary" tensor field on M defines a semibasic tensor field by the rule 
T T o 7r, so that the space of tensor fields on M can be treated as embedded in 
the space of semibasic tensor fields. 
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For a semibasic tensor field (T^ 1 '*j)(a;, £), the horizontal derivative is defined by 
Q . . d 

j!...j B \k ~ Q x k jl — 3s ~ kq^ Q£p 



ip jl-i _ 2^ fern jl-im-lMm+1- 



m — 1 

the vertical derivative by 



■ d ■ 

rpX\...% r rriZi...Z r 

jl...j s -k ~ g£fc ji...j 3 



and the modified horizontal derivative by 

rpi 1 ...i r _ rpi x ...i T _,_ \ f\yjrpii...i r 
± ji...j,:k - 1 j 1 ...j s \k \^\ I k 1 h-j s -j- 



The operators 

V, : C°°(#M) - C°°(0 r s+1 M), V. : C°°(#M) - C°°(^ +1 M), 

and 

V : :C~(/?rM)^C°°(^ +1 M), 

arc defined as 

(X7 , r TY 1 --- ir — T7 T»i~ ir ._ /jiti-ir 

and 

fY7 V7 rpi 1 ...i r . rpi 1 ...i r 

\ V -- 1 >h...j s k ~ y -- kl ji...] 3 — 1 ji...] s :k- 

For convenience, we also define V', V, and V : as 

In [301 133], the operators V| and V. were denoted by V and V respectively. 
Given u G C°°(TM \ {0}), we define 

Xu(x,0 = Cu:i = C^i+miu.j). 

Note that X restricted to SM coincides with G M . 
For V = (V*) G C°°(/3oAf), we set 

/t. m u 

chW:=V^ divV:=V3, divV :=Vl, (XVy = £ fe ^. 
Note if 7 is a magnetic geodesic, then 

(Xy)( 7 W,7(i))=V^(y( 7 (i),7(«))- 



Given a function it : SM — > IR, we will also denote by it its extension to a 
positively homogeneous function of degree on TM \ {0} (hoping that this will 
not yield any confusion). For a smooth <f> the smoothness properties of u defined 
by ()5.2p are determined by those of £(x,£,). As mentioned in the beginning of 
Subsection E21 the latter function is smooth on SM \ S(dM). All points of S(dM) 
are singular for £ as some derivatives of £ are unbounded in a neighborhood of such 
a point. Nonetheless, some derivatives are bounded. Let p be a smooth function on 
M such that dM = J o _1 (0) and | gradp = 1 in some neighborhood of dM. Define 

VpU(a:,£) = V : u(a;,£) - (V : u(a;, £), gradp(a;)) gradp(x). 
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Note that V^u(x,£) is completely determined by the restriction of u to the level 
set of p that contains (x, £). 

Lemma 5.1 (cf. [341 Lemma 3.2.3]). The semibasic vector fields Vil and V£ are 
bounded on SM\S{dM). 

Proof. Clearly, \7' p p = 0. Let h(x, £,t) = p(7 x ^(t)). Since /i(a;, £, 0) = p(x), it 
follows that V' p h(x, £, 0) = and therefore 

Vph(x,£,*) = a{x^,t)t 

for some smooth field a. Differentiating the equality /i(a:, £, £(a:, £)) = for (x, £) 6 
SM\S(dM), we obtain 

V- /i(x,^(z,0) + ^(x,^^,0)V^(x,0 

= V^(x,e,f(x,0) + (gi'adp(zj),77)V p ^,0 = 0, 
where (y,rf) = ijj^ x &(x,£) e d-SM. Hence, 

V A X & = -, ink \ 

(gradp(y),7?) 

Since £(x, £) < |L(y, 7/)|, boundedness of V p £(x, £) now follows from Lemma \2. 31 
Boundedness of V'£ is established similarly. □ 

5.2. Identities. The next identities are particular cases of those in [12j Lemma 4.6 
and Lemma 4.7]. For (x,£) <E SM: 

m 

2(V : u, V'(Xtt)) = | V : u| 2 + X((V : m, V u}) - div ((Xti)V'u) 

+ div ((X«)V : m) - (R £ (Vu), Vm) + (y(V'a), V : u), (5.3) 
2(V : u, V(Xu)> = |V : u| 2 + |V(Xu)| 2 + <Y(f), Vm) 2 - |X(Vu)| 2 , (5.4) 

where R ? (Z) = (i?j fcI f£'Z fe ) , 

= *j* + - Yi lk ) + gjMW - Y£Y t % 

and {R l jki) is the Riemann curvature tensor. 

For a unit speed magnetic geodesic 7, the operator C on smooth vector fields 
along 7 is defined by (see rT2J (45)]) 

C(Z) = R(j,Z)^-Y(Z)-(V z Ym, 

where Z stands for the covariant derivative along 7, Z = V-yZ. 

The operator C on smooth semibasic vector fields on TM \ {0} is defined by (see 
pi Subsection 4.6]) 

c(z)(x, £) = r*(z) - ypcz) - (v, z y)(0, 

where Rg(Z) = (^R l jki& £ l Z k ^J is the curvature operator. 
Notice that 

C(Z)( 7 (i),7(«))=C(Z(7(t),7(t))). 

In view of this identity, we henceforth omit the tilde in the notation of C, hoping 
that no confusion will arise. 
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Lemma 5.2. The following hold for (x,£) £ SM: 

|X(Vu)| 2 - (C(V'«),Vu> - (y(e),Vu> 2 = |V(Xu)| 2 -div((Xu)V : «) 

- div ((V : m, Vu)£ - (Xu)Vti) + div ((V : u, V'u}F(0), (5.5) 

|X(Vu)| 2 - 2(C(V'u), Vu) - (Y(£), V«) 2 

= |V(Xu)| 2 - 2div((Xu)V : u) - |V : u| 2 + 2(V : u,V(Xu)) 

- 2div ((V : u, Vu)£ - (Xu)Vti) + div ((V : u, Vm)7(0). (5.6) 

Proof. Note that 

X((V : «,Vu))(i,0 =e«V : n,Vu)) :j = ((V : ii,Vt 1 )f') :j - (V ! u,Vu)$ 

= div ((V : u, V«)f) - (V : m, Vh)Y/ = div ((V : it, Vu)£), (5.7) 

because Y? = by the skew-symmetry of Y. 
Using (|5.7[) , we first change (|5.3f) to 

2(V : u, V(Xu)> = |V : u| 2 +divV + div((Xu)V : u) 

- (Rf(Vu), Vu) + (Y(V'u), V : u), (5.8) 

with 

V= (V : u,Vn)£- (Xu)Vu. 

Next, 

div V = = Vj- + Y/KJ = div V + div (Y(V)) 

= divF + div((V : u, V'u)y(O) -div((Xu)Y(Vu)) 

and 

div((Xu)Y(Vu)) - ((X^Y/W^ = {Xu). J Y?g ik u. k + (Xu)g*Y>u. k . j 
= (V(Xu),Y(Vu)), 

because g kl Y/u.kj — by the skew-symmetry of Y and symmetry of mixed deriva- 
tives of u. 

Then flUHt takes the form 

2(V : u, V(Xu)) = |V : u| 2 + divU + div((V : u, V'u)Y(£)) - (V'(Xu), F(Vu)) 

+ div ((Xm)V : u) - (R c (Vu), Vu) + (F(Vti),V : u). (5.9) 
We have (see [HI Subsection 4.6]): 

(C(V«),Vu) = (R € (Vu),V'u) + (V'(Xu),y(V'u)> - (V : u,F(Vb)) 

+ <y(0,v u } 2 -<(V| (v . u) y)(0,v u ). 
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Also, 

(jz(v«,oe.v-«) + (v f y(v«), v«> - (v ( v„)no,v«) 
+ (y(0,0(^(vn),v u ) - (y(vu),o<y(0,v«) 
(R 5 (V ■«), v •«) - ((v (v . u) r)(e), v«> + (y(0, v u) 2 . 

Therefore, 

(R f (Vtt),Vu) = (C(Vu),Vu) - (v(Xu),y(Vu)> + (v : M,y(Vu)). (5.10) 

Using (|5TTU|) in f[579]) . we get 

2(V : u, V(Xu)> = |V : u| 2 + dWV + div((V : u, Vu)Y(0) + div ((Xu)V : it) 

- (C(V'u), V'u). (5.11) 

Subtracting (|5 . from (15. 4|l . we arrive at (I5.5[) . Multiplying (|5.11[) by 2 and 
subtracting the result from (15. 4[) , we arrive at (|5.6|) . □ 

5.3. Linear problem for 1-tensors. 

Theorem 5.3. Let (M, g, a) be a simple magnetic system, v a square integrable 
l-form, and tp a square integrable function on M . If the magnetic ray transform 
of the function <p(x, £) = Vi(x)£' + ip(x) vanishes, then ip = and v = dh for some 
function h G Hq(M). 

Proof. Recall that the space C 2 (M) consists of the pairs [v, <p], where v is a square 
integrable l-form and ip is a square integrable function on M, furnished with the 
norm (|3.16|) (see Section |3~3|) . In C 2 (M), we consider the subspace SC 2 (M) of 
solenoidal pairs [v, <p\, defined by the condition 5v = 0, and the subspace VC 2 (M) 
of potential pairs [dh, 0], where h G Hq(M). Then C 2 (M) decomposes as the 
orthogonal direct sum of SC 2 {M) and VC 2 (M). 

Associating each pair [v,(p] € £ 2 (M) with the function <j>(x,£) = Vi(x)^ +ip(x), 
define 

I[v,<p] = 1(f). 

Clearly, X vanishes on VC 2 (M) and, to prove the theorem, it suffices to show 
that kerJ nS£ 2 (M) = 0. 

By considerations similar to those in Section 4, one can prove, as in part (a) of 
TheoremSSl that KerX n SC 2 {M) C C l (M). 

Thus, we may assume under assumptions of the theorem that [v, <p\ € C l (M). 

Define u : SM — > M by means of (|5.2p . Then u satisfies the boundary value 
problem (|5.ip . As before, we preserve the notation u for the extension of it to a 
positively homogeneous function of degree on TM\{0}. Then for (x, £) G TM\{0} 
we have 

Xu(x,0=V i (x)e + \Z\l P (x). 

Now, we wish to integrate identity (|5.5[) over SM. However, u has singularities 
on T(dM), and we need some precautions against them. We proceed in the same 
way as in the proof of Theorem 3.4.3 in [34]. Let p : M — > R be a nonnegative 
smooth function such that dM = p _1 (0) and | gradpj = 1 near dM. For e > 0, let 



(R £ (Vm),Vh) = 
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M s = {x e M : p(x) > e}. Then u is smooth on SM £ . Integrating (|5.5p over SM e , 
we get 

{|x(v u)\ 2 - (c(v u), v •«> - (y(0, v u} 2 } ds 2 "- 1 ^^) 

= / \V (Xu)\ 2 di: 271 - 1 - [ div((Xw)V : M)dE 2 "- 1 

JSM e JSA/ e 

-/ div((V : u, V u)£ - (XM)Vu)dS 2 ™- 1 

■/SA^ 

+ / div((V : u, V'^K))^ 2 "- 1 ^^). (5.12) 

J SM e 

Now, we transform the last three integrals on the right-hand side by using the 
Gauss-Ostrogradskii formulas of f5U Theorem 2.7.1]. They give: 

div((Xw)V : M)dS 2 "- 1 = n t ({Xu)V : u,£_)dJ: 2n - 1 (x,0 

SMe JSM E 



n 



f (X«) 2 dY? n - ' 

JSM, 



because (V : u,£) = X«, 

div «V : u, Vu)( - (Xu)Vti) 

SM C 



2n-l 



/ {{V-u,Vu)Z - (Xu)Vu,gradp(x)} dZ^-'ix^) 

Jd(SM E ) 



(-1)"- 1 f ((V>,V«)(^gradp> - (V^,0<V-«,gradp» dS 2 



->2n-2 

because a straightforward calculation gives 

((V : it, V'it)C-(Xu)V'it,grad/9(a;)) = (V p u, V'u) (£, gradp) - (V p u, f)(V u, gradp), 
and 

div ((v-u, v«>y(f)) ds 2 ™- 1 ^, 

= (n-l) / («v i u,v«>y(0,0ds 2n-1 (*.0 = o 



because (y (£),£) = 0. 

Hence, f|5 . 1 2[) takes the form 



{|X(Vu)| 2 - (C(Vtt),Vu) - (y(£), V u} 2 } dE^-^as.O 



f \V (Xu)\ 2 cK?"- 1 -n [ {XufdY? 

JSM, JSM, 



,,v„; 2 ^ 2 ™- 1 

SM e «/SM e 

2n-2 

3(SM e 



-1)" / «V>,Vu}(£,gradp) - (V^,0(V-u,gradp)) dS 2 "" 2 (5.13) 
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Observe that by Lemma I5TT1 and (|5.2|) . the derivatives V p u and V'lt are bounded 
in SM \ S(dM). This fact allows us to pass to the limit in (|5.13p , arriving at the 
identity 

/ {|X(V u)| 2 - (C(V u), V u) - (Y(0, v u) 2 } <n? n ~\ x ,o 

J SM 

= [ \V (Xu)\ 2 di: 2 "- 1 - n [ (Xm) 2 ^ 2 "- 1 

JSM J SM 



+ (-1)" / ((Vj,u,Vu><e )g radp) - <V pU ,0(Vu,gradp)) dE 2 ™ 2 

d(SM) 



Note that the last integral on the right-hand side vanishes, because u\q(sm) = 
and hence V'u\q(sm) = 0. 

For € SM, Xu(x,£) = Vi(x)C + ip(x). As follows from Lemma 4.5.3 of 

[53] (and is easy to check directly), 

\V( Vl (x)C + <p(x))\ 2 dS 2 ™- 1 < n [ (Vi(x)C + ip(x)) 2 dE 2 ™" 1 . 

SA/ J SM 

Setting Z — V 'u, we thus have 

/ {|XZ| 2 - (C(Z), Z) - (Y(£), Z} 2 } dE 2 ™" 1 < 0. 

./SM 

By Santalo's formula (|A.4j) and the Index Lemma IA. 101 of Appendix A, this is 
possible if and only if Z = 0. So 1/(2, £) is independent of £, which clearly implies 
the sought conclusion. □ 

5.4. Linear problem for 2-tensors. For a magnetic system (M, g, a) and (x, £) G 

■SM, put 

fc^z, = sup{2K(x, a s>n ) + (Y( V ), £} 2 + (n + 3)|y(?7)| 2 - 2((V^)(0, r?}}, 

where the supremum is taken over all unit vectors 77 £ T^M orthogonal to £, and 
K{x,a^^) is the sectional curvature of the 2-plane cr£ i?) spanned by £ and 77. 
Define 

= max{0,/c(x,^)} 

and 

k(M,g,a) = su P T 7 / 7 fc+( 7 (i) )7 (f)) dt, 



where the supremum is taken over all unit speed magnetic geodesies 7 : [0, T 7 ] — > M 
running between boundary points. 

Theorem 5.4. If {M,g,a) be a simple magnetic system with k(M,g,a) < A, then 
I is s-injective. 

Remark 5.5. Note that the condition k(M,g,a) < 4 holds if (M, g) is negatively 
curved and the C 1 -norm of Y is small enough. Also, it is easy to see that this 
conditions is valid for every sufficiently small simple piece of any magnetic system. 

Proof. Assume that [h,/3] £ ker/ n «SL 2 (M). Then, by Theorem gTBJa) , the pair 
[h, (3] is smooth. Define u : SM — > R by means of (|5.2[) . Then u satisfies the 
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boundary value problem (|5.1j) . As before, we preserve the notation u for the exten- 
sion of u to a positively homogeneous function of degree on TM \ {0}. Then for 
(x, £ TM \ {0} we have 

Xu = £ '/^(^ee'+^We- (5-14) 

Now, we wish to integrate identity (|5.6[) over SM. Again u has singularities on 
T(dM). We overcome this obstacle in the same way as in the proof of the previous 
theorem. Thus, we deduce 

f {\X(Vu)\ 2 -2(C(Vu),Vu)-(Y(0,S7u) 2 } d^-^x^) 

J SM 

= [ \v {Xu)\ 2 a?? 11 - 1 -2 f div((Xw)v : u)ds 2 "- 1 - f iv^i 2 ^ 2 "- 1 

JSM JSM JSM 

+ 2 [ (V-u, V{Xu))dY, 2n - 1 - 2 / div((V : u,Vu}£- (Xw)Vu) dS 2 ™" 1 

J SM J SM 

r i) 

2n-l 



+ / div((V : u, Vu)Y(£))dT,' 

JSM 

Again we use the Gauss-Ostrogradskii formulas to transform integrals of diver- 
gent form. Denoting the leftmost side of the above formula by A, we obtain 



A= [ \V {Xu)\ 2 dY, 2n - 1 - j {2n(Xu) 2 + |V : u| 2 } d£ 

JSM JSM 

+ 2 [ (V^VtXu)}^ 2 "- 1 

JSM 



SM 

<[ \V(Xu)\ 2 dZ 2n - 1 - (2?i + l) f {XufdY? 11 - 1 

JSM JSM 

+ 2 [ (V : u,V{Xu))dT l 2n -\ (5.15) 

JSM 



'SM 

where we have used the inequality \Xu\ = |{V : u,£)| < |V : it 
From (|5.14j) we have for G SM 

(Xu) 4 = 2h lj e -Zlhije? 

Therefore, 



/ \V {Xu)\ 2 dY? 11 - 1 = [ \2h(0 - (h,^ 2 )^ + (3\ 2 dY 2 

JSM JSM 



-"«.) ■ " s"-<. ; < ->■ ' 1 ' A '"" 1 

ism Jsm 

{A\h(i)\ 2 - 3(h,e) 2 + \(3\ 2 +4(h(0,P) ~ (h,e)(P,0} dZ 2 "- 1 



SM 

{4\h(0\ 2 ~Hh,e) 2 + \(3\ 2 }dJ: 2n -\ (5.16) 



SM 



f {XufdY? 11 - 1 = f {(/ l ,e 2 ) 2 + (/3,0 2 }^ 2 "- 1 . (5.17) 

JSM JSM 

Now, we transform the last integral in (|5.15p . Once 

V : it = Viit - Y(V.u), 



it follows that 

(V : u,V(Xm)) = (Vi«,V.(Xu)) - (Y(V.u), V.(Xu)}. (5.18) 
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Observe that 

(V|u, V.(X«)> = g kl u\ k {2hi :j e ~ tihiit'e + Pi) 
= 2{ug M h l3 ^ k - 2ug kl h lj>k ? - (u^*/iy^')| fc + f ) ]k + (u/3% - u(3 k k 

= div (W) - 2u(5h) j £ j + uG{h tl Ci 3 ) - u6{(3) 

= d\v(W) - (n-l)(uY((3),0+uG((h,e)), (5.19) 

with 

W k = 2ug kt hij^ - <%rr + u(3 k . 

Here we have used the equations p.23[) for solenoidal pairs. 

Integrating (|5 . 19[) and transforming integrals by means of Gauss-Ostrogradskh 
formulas, we get 



f (V|7!,V.(Xli)}dI] 2 "- ' 
J SM 



-(n-l) (uY(/3),0dX 2n - 1 + uG((h,Z 2 ))dZ 2n - 

JSM JSM 
SM J SM 



(Y(0),Vu)dY; 2n - 1 - G M (u)(h,£ 2 )dE 



2n-l 

SAf J SM 

+ / (rr/^-)^^ 2 )^ 2 "- 1 



SM 



(FfVu),/?)^" 1 - / (/i,£ 2 ) 2 rfS 2 ™" 1 

SAf JSA/ 



(V Ul F(0)(^r}^"-\ (5.20) 

SM 



Next, 



(F^.^.V.CXu))^ 2 "- 1 

9M 

= / (Y(V.u), (2hije ~ ZihrjCe + A)} dY? n - 1 

JSM 

= 2 ( (Y(Vu),h(0)dZ 2n - 1 + f (Vu,Y(0)(h,e)dZ 2 

JSM JSM 



+ / (Y(Wu),P)dE 2n - 1 , (5.21) 

JSM 



where h(£) = {</'/',,:' i. 
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From IpUS]) . (j5~20l) . and dOTI) we get 



f <V : u,V(Xw))d£ 2 

JSM 



->2n-l 

\ V 11, V ^YV'U M CI/ 

'SM 

t2\2 jv2n-l 



(/i,r> 2 dS 2n_1 -2 / (r(Vu),/i(0)dS 2 ™- 
<-/ (h,^) 2 dE 2n - 1 + ^— [ \h(0 1 2 dE 2 "- 1 



SM n + 2 J SM 

n + 2 



/ lytVajprfE 2 "- 1 . (5.22) 



ISM 

Using (j5TT6l) . (|5T7| . and (IE22| in (J53U), we find that 
A< f {4\h(£)\ 2 -3(h,Z 2 ) 2 + \f3\ 2 }dZ 2n - 1 

J SM 

-(2n + l) / {(/i,e 2 ) 2 + (/3,0 2 }^ 2 "- 1 -2 / (^e 2 ) 2 ^ 2 "" 1 

JSM JSM 

2 jv^2n— 1 i / i r,\ I IV/T7'».M2 jv2ti — 1 



n + 2 
4(n + 3) 

n + 2 JSM 



|/ l (0rrfS^- i + (n + 2) / |r(Vu)|"dS 2 

SM JSM 



^(OpdS 2 "- 1 -2(n + 3) / {h,?) 2 dY? 11 - 1 

SM 



+ / l/?! 2 ^ 2 ™" 1 - (2n + l) / {13,0 d^ 

JSM JSM 



+ (n + 2) / \Y(Vu)\ 2 dT, 

JSM 

As follows from Lemma 4.5.3 of [53] (and is easy to check directly), 

WOfdE 2 "- 1 < / (/i.C 2 ) 2 ^ 2 "- 1 , 

2 y SM 

2 dY 2n - 1 <nj (/3,0 2 d^ 2n -\ 

SM JSM 

Therefore, 

|2 jv^2n-l 



2 j^2n-l 



A < (n + 2) / |F(Vu)| 2 dS 2 

JSA/ 

Setting Z = V'u, we rewrite this as follows: 

{|X(Z)| 2 - 2{C(Z),Z) - (Y{Z),0 2 - (n + 2)\Y(Z)\ 2 } dS 2 ™" 1 < 0. (5.23) 

ISM 

We need the following lemma whose proof is given below. 
Lemma 5.6. Let 7 : [0,T] — ► M be a unit speed magnetic geodesic. If 

[ T k+( 7 (t)^(t))dt<A/T, 
Jo 

then for every smooth vector field Z along 7 vanishing at the endpoints of 7 and 
orthogonal to 7 we have 

jT {\Z\ 2 - 2(C(Z),Z) - (Y(Z),j) 2 - (n + 2)\Y(Z)\ 2 } dt > (5.24) 
with equality if and only if Z = 0. 
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Continuing the proof of the theorem, observe that, by Santalo's formula (|A.4|) 
and Lemma [5.61 inequality (|5.23j) may hold if and only if Z = 0. So Vu = 0, 
which means that u(x,£) is independent of £. Since the pair [h, 0\ is solenoidal, this 
readily yields the sought conclusion. □ 

Proof of Lemma [57h\ Denote the left-hand side of (I5.24|) by B. Since 
C(Z) = R^Z)-Y(Z)-(V z Y)(j), 

we have 

f-T 



D 



[\Z\ 2 - 2(R^(Z), Z) + 2(Y(Z), Z) + 2<(V Z Y)( 7 ), Z) 

(Y(Z) } j) 2 -(n + 2)\Y(Z)\ 2 } dt. 



Let e,, i = 1, . . . , n be an orthonormal frame at 7(0), with e n = 7(0). Let Ei(t), 
< t < T, be a vector field along 7 satisfying the equation 

^ - Y{Ei) 

with the initial condition Ei(0) = e^. Surely, E„ = 7. Since 

|(^,^i> = (Bf,^-) + - {Y{E t ),Ei) + {Ei,Y(Ej)) = 0, 

we see that Ei(t), . . . , is an orthonormal frame for each t. 

Consider the expansion of Z in this frame: 

Z(t) = z K {t)E K (t), «6{l,...,n-l}. 

Then 

Z = z K E K + z K E K = z K E K + z K Y{E K ) = z K E K + Y{Z), 
\Z\ 2 = J2 z K z K + 2Y kX z k z x + \Y(Z)\ 2 , 

(Y(Z), Z) = (z K Y(E K ) + Y 2 (Z),z x E x ) = -Y kX z*z x - \Y{Z)\ 2 

with 

Y kX = (Y{E K ),E X ). 

Therefore, 

B = I { £ z K i K [2K{a^ z )\Z\ 2 + {Y(Z)^) 2 

+ (n + 3)\Y(Z)\ 2 -2{(\7 z Y)(i),Z)]}dt 
> {£i K i K -k ll \Z\ 2 }dt = £{i K i K -fc A1 z K z K }dt. 
Now, the claim follows from Lyapunov's inequality [16[ p. 346]. □ 
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6. General rigidity theorems 

6.1. Rigidity in a given conformal class. We first state a rigidity theorem, fix- 
ing the conformal class of a metric. The theorem below generalizes the correspond- 
ing well-known result for the ordinary boundary rigidity problem, see [9j 1231 124] . 

Theorem 6.1. Let (g,a) and (g',a') be simple magnetic systems on M whose 
boundary action functions A\qmx8ai and A'\dMxdM coincide. If g' is conformal 
to g, i.e., g' = uj 2 {x)g for a smooth positive function uj on M, then u) = 1 and 
a' = a + dh for some smooth function h on M vanishing on dM , hence (</, a') is 
gauge equivalent to (g,a). 



Proof. In view of Lemma 12.11 ui = 1 on the boundary of M. Next, using Lemma 
we see that the scattering relations S and S' of both magnetic systems coincide. 
Let us show that ui = 1 on the whole of M. Note that for (x, £) £ d+SM 

A'« c ) < A'( 7 *,s) = \j ( ' ?) ^ 2 ( 7 ,,« (*)) dt + ht(x, 0-1 (6-1) 



7*,e 



Using (|2.9|) and (|2.8|) , whence we obtain 

— / 

Wn-x Jd + SM 



V oV(M) = ^-/ A'( 7 ^)d M ( s ,0 



1 



(6.2) 



^n-l Jd+SM Uo J 

= i^ W 2 dVol+ivol 9 (M). 
On the one hand, by Holder's inequality 

cj 2 dVo\<\[ uj n dVo\\" \ [ dVol\ = VoV(M)» Vol s (M)™ 
M 17m J 17m J 

with equality if and only if u> = 1 . 
It follows that 

VoLy(M) < 1 VoV(Af)^ Voi g (Af)^ + i Vol fl (M). (6.3) 

However, by Theorem ES VoL/(M) = Vol a (M), which implies that ^673j) holds 
with the equality sign. This means that (|6.2p holds with the equality sign. Thus, 

LU=l. 

Now, (|6.ip and the equality A' (7^) = A(j x ^) yield 

a > / a'. 

In view of (|A.4[) and (|2.8p . we conclude now that 

(a' - a) = 

for all (x, £) £ d+SM. By Theorem [531 we see that a' — a = dh for some function 
h on M vanishing on dM. This completes the proof of the theorem. □ 
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6.2. Rigidity of analytic systems. 

Theorem 6.2. If M is a real-analytic compact manifold with boundary, and (g,a) 
and (g' , a') are simple real- analytic magnetic systems on M with the same boundary 
action function, then these systems are gauge equivalent. 

Proof. Consider the diffeomorphism / : M — > M constructed in the proof of 
Theorem 12.21 Recall that in a neighborhood of dM, f coincides with the map 
exp aM o (expg M ) _1 , which is obviously analytic for analytic systems. 

As shown in that proof, g and f*g' have the same jets on dM. Hence, by an- 
alyticity, g and f*g' coincide in some connected neighborhood of dM, and / is 
thus an analytic isometry of a neighborhood of dM of the the analytic Riemannian 
manifold (M, g) onto a neighborhood of dM of the analytic Riemannian manifold 
(M,g'), which is moreover the identity when restricted to dM. Now, the same 
arguments as in the proof of [HJ Theorem C(a)] show that / extends from a neigh- 
borhood of dM to an analytic isometry / of (M,g) onto (M,g'). Now, the sought 
conclusion follows from Theorem 16 . 1 1 applied to the magnetic systems (M,g, a) and 



6.3. Rigidity of reversible systems. A magnetic system (M,g,a) is said to be 
reversible if the flip (x, £) i— > (x, — £) conjugates tp* with tp^ 1 . It is easy to see that 
a system is reversible if and only if da = 0, i.e., if and only if the Lorentz force Y of 
the system vanishes. In this case magnetic geodesies are nothing but the ordinary 
geodesies of the Riemannian manifold (M,g); moreover, simplicity of a reversible 
magnetic system (M, g, a) is equivalent to simplicity of the Riemannian manifold 



It is interesting to know whether reversibility of a magnetic system can be es- 
tablished by boundary measurements. Observe that if a system is reversible then 



for all (x,£) G d + SM. Let us call the systems satisfying (|6.4I) boundary reversible. 

Theorem 6.3. A simple magnetic system is boundary reversible if and only if it 
is reversible. 

Proof. Let (M, g, a) be a boundary reversible simple magnetic system. Consider 
the magnetic system (g, —a) on M. Note that if 7 : [0, T] — > M is a unit speed 
magnetic geodesic of the system {g,a), then the curve 7 : [0,T] — ► M, defined as 
7(i) = j(T — t), is a unit speed magnetic geodesic of the system (g, —a). 

For x, y £ M, let j x<y : [0, T x , v ] — > M denote the unit speed magnetic geodesic 
of the system (5, a) from x to y. By Lemma I A. 51 



(M,/y, /*«')• 



□ 



(M, 5 ). 



<S(-S(x,£)) = (z,-0 



(6.4) 



A(7x, a ) < A(7j / , :e ), 



i.e. 




(6.5) 



Interchanging x and y, we receive 




(6.6) 
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From (|6.5|) and (|6.6|) we get 

a + I a>0. (6.7) 



'7 H ,x 

For (x,£) 6 d + SM and y = r y x ^(£(x, £)), we have in view of (16. 4p 

7x,s/=7x,e. 7i/,a:=7s,£, 

where 7^ : [0,£(x,£)] — > M is the unit speed magnetic geodesic of the mag- 
netic system (M,g, —a) with initial conditions 7 x ,f(0) = x, 7x^(0) = £, and with 
7x,£W>,£)) = y- Therefore, (JBT7]) yields 



a- a>0. (6.8) 
By Santalo's formula (j A.4|) and (|2.8[) . we have 

/ W a ~ / 4 ^,0=0. (6.9) 

Jd+SM [J-y^ J 

Combining (|6.9[) and (|6.8p yields 

/ a = / a (6.10) 
for all (x,£) G d+SM. From (16. 5|) we then obtain 

t( x ,z)<i{x,s), 

and from (16.61). 



Thus, 



i(x,z)<i{x,o. 

i(x,Q=t{x,t). (6.11) 



Now, by ([00]) and (j6TTj) 

A(7c,s) = Mlx.z) = A(x,y). 

This implies that is a unit speed magnetic geodesic of the system (g, a) which 
joins x to y. By simplicity, we then have 

%,S = 7*,«' 

Valid for all (x, £) e d+SM, this equality means clearly that the magnetic system 
(M, g, a) is reversible. □ 

Theorem 6.4. A simple reversible magnetic system (M, g, a) is magnetic boundary 
rigid if and only if the simple Riemannian manifold (M, g) is boundary rigid. 

Proof. Suppose that a simple reversible magnetic system (g, a) on M is magnetic 
boundary rigid. Let g' be a simple Riemannian metric on M whose boundary 
distance function equals the boundary distance function of g. We must prove that 
g' is isometric to g by an isometry which is the identity on the boundary. 

Without loss of generality we may assume (see, e.g., [2Ql Theorem 2.1]) that 

g'\dM = g\dM- (6-12) 
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Then equality of the boundary distance functions of g and g' and (|6.12j) are well 
known to imply equality of the scattering relations of the Riemannian manifolds 
(M,g)axid{M,g') (see, e.g., M)- 

Consider the magnetic system (</ , a) on M . Since da = 0, this system is re- 
versible as well. It is trivial that for a reversible system, the scattering relation of 
the magnetic system coincides with the scattering relation of the underlying Rie- 
mannian metric. Therefore, the magnetic systems (<?, a) and (</, a) on M have 
the same scattering relations. In view of (|6.12[) and Lemma [231 the boundary ac- 
tions functions of these systems coincide. From the magnetic boundary rigidity of 
(M, g, a) we infer that the metric g' is isometric to g by an isometry which is the 
identity on the boundary, as required. 

Now, let (M, g, a) be a simple reversible magnetic system such that the simple 
Riemannian manifold (M, <?) is boundary rigid. Let (g',a') be a simple magnetic 
system on M whose boundary action function A'\gMxdM equals the boundary ac- 
tion function A\gMxdM of the magnetic system (g,a). By Theorem I2.2i we may 
assume that g'\dM = g\dM and a'\gM = ol\qm- Now, Lemma 12.51 implies that the 
scattering relations of the magnetic systems (M, g, a) and (M, g' , a') are the same. 
Since (M, g, a) is reversible and hence boundary reversible, (M, g' , a') too is bound- 
ary reversible and hence, by Theorem l6.31 it is also reversible. This implies readily 
that the scattering relations of the Riemannian manifolds (M,g) and (M, <?') are 
the same. Applying Lemma [2.61 to the magnetic systems (M,g,0) and (M, g',0), 
we see that the boundary distance functions of the metrics g and g' are the same. 
Hence, by the boundary rigidity of the Riemannian manifold (M, g), there is a dif- 
feomorphism / : M — > M, f\dM = identity, such that g' — f*g. Then Theorem 
16.11 applied to the magnetic systems (M, g, a) and (M, f*g', f*a'), gives the gauge 
equivalence of the magnetic systems (g,a) and (g',a'), as required. □ 

6.4. Generic local boundary rigidity. We will prove that near each (go,ao) in 
the generic set Q k of Definition 14. 101 the action on the boundary determines (<?, a). 
Note that, by Theorem 15.41 this generic set Q k contains the magnetic systems 
(g,ct) with k(M,g,a) < 4, in particular magnetic systems in which the underlying 
Riemannian metric is negatively curved and the magnetic field is sufficiently small. 

Theorem 6.5. Let kg be as in Theorem \4-H\ There exists k > ko such that for 
every (go,cto) € Q k , there is e > such that for any two magnetic systems (g,a), 
(<?' , a') with 

11.9 ~ 5o||c fc (M) + ||a - ao||c fc (M) < e, \W - 9a\\c k (M) + \\ a ' - ^o\\c k (M) < £ 
we have the following: 

A s , a = Ag' Q ' on dM x dM (6.13) 

implies that (g,ct) and (g',a') are gauge equivalent, i.e., g' — ip*9> a ' — ip*ot + dcf> 
with some C k+1 (M)-diffeomorphism ip : M — » M, fixing the boundary, and some 
C k+1 function (j> vanishing on dM . 

Observe that g is solenoidal (with respect to itself), and if we replace a by its 
solenoidal projection a s :— a — d<j>, where A g 4> = 8a, 4>\dM — 0, then Sa s = 
as well. Therefore, [^g, —a s ] is not necessarily a solenoidal pair in the sense of 
Definition 13.71 instead 

1, 



"1 








-a s 




2 9 







-Y{a s ),0 



(6.14) 
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We will prove an analogue of [ITJ Theorem 2.1]. 



Lemma 6.6. Let (g,a) S C fc,A1 , < [X < 1, k > 2, be a simple magnetic system on 
M with 5a = 0. Then for any other (g',a') close enough to (g,a), there exists a 
magnetic system (g' ,a') gauge equivalent to (g',a') and satisfying (|6.14p , i.e., 







}~9',-a'_ 





-Y{a),0 



(6.15) 



where 8 and Y are related to {g,a). 

Moreover, if \\{g' , a') — (g, a)||c?fc,M < e with e -C 1, then \\(g', a') — (g, a)\\ci°,» < 
Ei, with £i = Ei(e) — > as e — > 0. 

Proof. Our argument is much the same as that in the proof of [111 Theorem 2.1]. 
So we merely sketch it. 

We have to show that there exist a diffeomorphism / of M, fixing dM, and a 
function tp, vanishing on dM, such that <?' = f*g', a' — f*a' + dtp, and 15[) holds. 
If / is close enough to the identity, we can identify / with a certain vector field v 
as follows. If v is a vector field vanishing on dM and |Vt>| 5 < 1, the map 

e v (x) = exp x (v(x)) 

is well defined on M, with image in M again, and it is a diffeomorphism for Hvllc^.c 
small enough. Then v i— » e v has an inverse defined as follows. If / is a diffeomor- 
phism close enough to the identity, set Vf(x) = exp~ 1 (/(a;)), i.e, Vf(x) = 7(0), 
where 7 : [0, 1] — > M is the geodesic such that 7(0) = x, 7(1) = f(x). The exis- 
tence of such a geodesic follows from the strict convexity of the boundary, following 
in turn from the simplicity assumption. Clearly, the differential of the map v 1— > e v 
at v = is the identity transformation. 

Let h = i(</ — g), [3 = —(a' — a). Condition (|6 . 1 5|) then takes the form 



7? — 1 71 — 1 

6 (e* v (g/2 + h)) + —Y ((e*(a + 0) + dip) = —Y{a), 

-S(e* v (a + (3) + dip) = 0. 
We define the map (recall that {g, a) is fixed) 
F([v,p],[h,f3}) 

— -Y ((e* v (a + !3)+ dp) , -5 (e*(a + f3) + dip) 



(6.16) 
(6.17) 



S(e* v (g/2 + h)) + - ^ 
Then (|6.16p , (|6 . 1 7[) can be rewritten as 
F([v,<p],[h,p}) = 



-Y(a),0 



(6.18) 



We want to solve this equation for [v, <p] if [h, 0] are small enough and g, a are fixed. 
If [h, 0} — [0,0], then [v,ip] = [0,0] is a solution by (|6.14[) . To show solvability for 
[h, 0] small enough, we apply an implicit function theorem in Banach spaces as in 
[IT] . For this purpose, we need to compute the derivative F' v ^([0,0], [0,0]). 
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Set [h, [3] = in IpTTg) and use (j3~20l) to deduce 
F( vM ([0,0},[0,0])[v,ip] 



Sd s v 



~-Y(Y(v) + d(v, a) + dip), -S(Y(v) + d(v, a) + dip) 



= Sd[v, -(v,a) - (p\, 

where v is considered as a 1-form by lowering the index. Thus, the derivative above 
is the superposition of the map [v, ip] ^— > [u, — (v, a) — <p], which is an isomorphism 
(and this map and its inverse preserve the zero boundary conditions), and Sd with 
Dirichlet boundary conditions, which is also invertible in appropriate spaces. We 
refer to [11] for more technical details. □ 

The next lemma states, loosely speaking, that gauge equivalent pairs differ by a 
potential one, modulo a quadratic term. 

Lemma 6.7. Let (g',a') and (g,a) be in C , k>2, and gauge equivalent, i.e., 

g' = ip*g, a' — ip*a + dip 

for some diffeomorphism ip fixing DM and some function ip be vanishing on dM . 
Set f = [5(3' — g), — (a' — a)]. Then there exists w, vanishing on dM , such that 

f = dw + f 2 , 

and for (g,a) belonging to any bounded set U in C k , there exists C{U) > such 
that 

\\h\\c--' < C{U)U - Idll^-x, ||w|| c *-i < C(U)U ~ Id|| c *-i. 

Proof. As in the proof above, set v{x) = ex.p~ 1 (ip(x)), which is a well defined vector 
field if ip is close enough to the identity in C 2 (it is enough to prove the claim in this 
case only) and v — on dM. Set i/j t (x) = exp x (rv(x)), < r < 1. Let g T — ip*g. 
Then the Taylor formula implies 



where 



1 



\h\ < — max 

2 r G [0,l] 



h = g + 2d s v + h, 



dT 2 



and 2d s v is the linearization of g T at t = 0. To estimate h, write 

dip*dip l T 

and differentiate twice w.r.t. r. Notice that 



8t 2 



<C\\v\\ 



8t 2 



<C\\v\\^. 



This yields the stated estimate for the first component of f~2 for k = 2. The estimates 
for k > 2 go along similar lines by expressing the remainder h in its Lagrange form 
and estimating the derivatives of h. 

The analysis of the second component is similar, using (|3.20p . In particular, we 
get that w = [v, —a(v) — <j)\ which corresponds well with the linearization formula 
(I3~20l) . □ 
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Proof of Theorem \6.5\ We can assume that (g' , a') is replaced by its gauge equiva- 
lent pair that satisfies the assumptions of Lemma l6.6l (where 8 and Y are related to 
(g,aj). Set f = [(g' — g)/2,—(a' — a)]. Then = 0. We assume that (g,a) belongs 
to a small enough neighborhood of (go,ao), so that (g,a) € Q k , and the constant 
C in Theorem 14.31 is uniform in that neighborhood. By the second statement of 
Lemma 16.61 and the assumptions of Theorem 16. 5[ 

||f||o*<ei(e), (6-19) 

where £\ — > as e — * 0, and k 3> 1 is fixed (this requires the original (g',a') to 
be in Q k+1 if we want to avoid the C k ^ spaces). Moreover, HgHc?* + ||o:||c? fe < A, 
where A > depends on (go, cxq) and on an upper bound £o of e. All constants C 
below will depend only on A and will be uniform in e < sq. 

Let ip be a diffeomorphism in M that maps boundary normal coordinates w.r.t. 
g into boundary normal coordinates w.r.t. g' near dM, then extended to the whole 
M. In other words, if i/'i that maps a neighborhood of dM into dM x [0,6] defines 
the scmigcodesic normal coordinates related to g and if ^2 is defined in the same 
way corresponding to g' , then tp — ^2 ipi' ■ Clearly, 

H - H||c*-> < C\W - gWc- < C\\f ||c*. (6.20) 

Set g' — ip* g' , a' — i\)*oi + dtp, where (p is such that in boundary normal coordinates 
54 = 0, see LemmaHU Set f = [g'/2, -a'] - [g/2, -a]. By Lemma[6J]and (f6T20j) . 

f = f + dw + h, ||h|| c *- a < C\\ff ck , (6.21) 

with w\qm = 0; so, roughly speaking, f and f differ, up to a quadratic term, only by 
a potential term. By Theorem 12.21 f vanishes on dM together with its derivatives 
up to any fixed order m, if k 3> 1. We will choose m below. We are going to prove 
that f = 0, and this would prove the theorem. 

Since A S . Q = A g r, a r — Ag^g' on dM x dM, by the linearization Lemma HOI 

||/f|U~ < cufii'x. 

By ([323, \\Nf[[ La o {Ml) < C||Zf|| £ -; therefore, 

||^f|U-(MO < CPIIc- (6-22) 

We apply Theorem 14.31 to get 

||f s |U*(M) <C\\Nf\\ mMl) , (6.23) 

where we estimated the H 2 norm by the H 2 norm, which is finite for f . Note that 
for the solenoidal projection f s of f we have, by (|6.21[) . 

P = f + h s ; (6.24) 

in particular, 

||f s |U*(M) > \\{\\lhm) - \\h s \\ LHM ) > ||f|U=(M) - C||f|| 2 c2 . (6.25) 

Using interpolation estimates in H s (Mi), see [5TJ Theorem 4.3.1/1], and (|6.22p . 

we get 

\\N~{\\ HHMl) < C||JVf||^ (Wi) ||JVf||^ < C'||f||^ 4/S , - > 2. (6.26) 
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To estimate the first factor in the middle term, we used the fact that d a i\dM = 
for \a\ < m and therefore, if m > s — 2, 

\\N~f\\ H . iMl ) < C\\f\\ H .- 1(M) < C'||f|| C3+1(M) < C" 

if k > s + 1. The last inequality follows by comparing f, f using (|6.21|) and the 
estimate on w in Lemma \Q.7\ combined with (|6.20|) . Combine ()6.23|) . (|6 . 25[> . and 
flHUj)) to get 

||f|| L2(M) <c(||f|| 2 cl 4/s + ||f|| 2 c; 

Using again the fact that ||f ||c fe (M) < C 1 (fc ^ 1), Sobolev embedding estimates, 
and interpolation estimates in H S (AI), one gets 

||f || L2(M) < c||f||i 2 -^ 

with < /i < 1 that can be chosen as close to 1 as needed, provided that k ^> 1. It 
is enough now to choose s > 2 and /x so that (2 — 4/s),u > 1 to deduce that f = 
if e < 1, see dfiHU). □ 



7. Rigidity of two-dimensional systems 

The main result of this section is the following rigidity theorem for two-dimensional 
systems: 

Theorem 7.1. If dimM = 2 and (g,a) and (g',a') are simple magnetic sys- 
tems on M with the same boundary action function, then these systems are gauge 
equivalent. 

This theorem generalizes the boundary rigidity theorem for simple Riemannian 
surfaces which was established in [35] . Our proof of Theorem 17.11 mimics that of 
the mentioned theorem in [32j . 

First of all, by Theorem 12.21 and Lemma 12.51 we may assume that g' coincides 
with g on dM, 

g'\dM = g\dM, 

and that the scattering relations of the magnetic systems (M, g, a) and (M, g'a!) 
coincide: 

S'=S. (7.1) 

The crucial step then consists in establishing that the scattering relation of a 
two-dimensional magnetic system (M, g, a) determines the Dirichlet-to-Neumann 
(DN) map associated to the Laplace-Beltrami operator of the Riemannian manifold 
(M,g). 

It is proved in [211 [19] f° r two-dimensional manifolds that the DN map determines 
the conformal class of the Riemannian metric up to an isometry that is the identity 
on the boundary. 

Afterwards, the proof of Theorem 17. II is finished by applying Theorem 16.11 which 
claims magnetic boundary rigidity within a given conformal class. 

Derivation of the connection between the scattering relation and the DN map 
is based on the properties of the magnetic ray transform and the commutation 
formula between the magnetic flow and the fiberwise Hilbert transform. 

We proceed with describing the needed properties of the magnetic ray transform. 
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7.1. More about the magnetic ray transform. We recall that, given a notation 
F for a function space (C k , L p , H k , etc.), we denote by J-(M) the corresponding 
space of pairs [v, tp], with v a 1-form and tp a function on M. Also, recall that in 
the space C 2 (M) we consider the norm (|3.16[) defined as 



Iwll 2 = 



(M 2 +<p 2 ) dVol. 

' M 

Associating each pair [v,<p] € C 2 (M) with the function <j>(x,£) = Vi(x)t? +<p(x), 
we may consider £ 2 (M) as a subspace of L 2 (SM). 

Consider the restriction X of the magnetic ray transform I to C 2 (M); 

£) = / <lx,dt))jl f (t)dt+ <p( lx>£ (t))dt 
Jo Jo 
= I lV (x, + I Q y{x, 0, (x, e 

By Lemma 13.21 we receive an operator 

1:C 2 {M) -^Ll(d+SM). 

Let 

J* : Ll(d+SM) -» £ 2 (M) 
be its dual. The same calculations as in Subsection 3.4 show that 

with 

\JS X M J JS X M 

The following holds: 

kerJ = 7'£ 2 (M), (7.2) 

Imrc5£ 2 (M). (7.3) 

Equality (|7.2p is just the claim of Theorem 15.31 Equality (|7.3[) follows, as soon as 
ImX* is in the orthogonal complement to kerl and the orthogonal complement to 
VC 2 (M) is SC 2 (M). 
Consider the operator 

Af = 1*2" : C 2 (M) — > C 2 (M), 
which can be written down as 

Af[v,(p] = [Afnv + J\f w ip,Af iv + AfoQtp}. 
Considerations similar to those in the proof of 14. II show the following: 
Proposition 7.2. M is a ^DO in M mt of order —1 with principal symbol 
op (AO = diag (op(A/ii), v p (Nm)), 

a p (Af 11 )i(x,0=c n \^\- l (Si-e^/\e), 
W /iere |£| 2 = 0«(s)&& and ^ = <^(x)&. 
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Aw(x, £) 



Define 

C™(d+SM) = {we C°°(d+SM) : w B 6 C°°(SM)}, 
where ufl is the function that is constant along the orbits of the magnetic flow and 
equals w on d + SM . This space can be described in terms of the scattering relation 
S alone; by the forthcoming Lemma ITUl 

C™(d+SM) — {w G C°°(d + SM) : Aw G C°°(d(SM))}, (7.4) 

where 

'w(x,0, {x,Qed+SM, 
v too5- 1 (x l f)i efl-SM. 
The following is an analog of [29l Theorems 3.3.3, 3.3.4] (see also [32, Theorem 
1.4], [31, Theorems 4.1, 4.2]). 

Theorem 7.3. Let (M,g,a) be a simple magnetic system. Then, for every pair 
[v, ^} G C°°(M), there exist w G C™(d + SM) and f G C*°°(M) such that 

[ v + V/, ip] = Tw. 

(Note that if v is solenoidal then, in view of (|T. 3[) . / is harmonic in M .) 

Proof. Our argument is the same as in [29l E21 El] . Embed M into a closed manifold 
M and extend g and a smoothly to a Riemannian metric g and a 1-form a on M. 

If U C M is an open neighborhood of M with smooth boundary, then (U,g,a) 
is also a simple magnetic system if dU is close enough to DM. Henceforth such an 
U is assumed to be fixed. 

Denote the magnetic ray transform for (U,g,a) by Ixj, and denote by tm the 
operator of restriction to M. We have the following analog of Theorems 3.3.1, 
3.3.2] (see also [32 Theorem 3.1], [31j Theorem 4.3]): 

Lemma 7.4. For every pair [v,ip] £ Ti. s+1 (M), s > 0, there exists a pair [u,ip] G 
H S (U) and f G H S+2 (U) such that 

[v + Vf, ip] = r M IuIu[u, i>\- 

Again, the proof is similar to the one in 29, 32, 31]; we will merely sketch it. 

Cover M by finitely many open sets Uk such that U = U%, D M = for 
k > 2, and (Uk,g, ct) is a simple magnetic system for every fc. Let {/i^} be a subor- 
dinate partition of unity such that h\\M — 1. Consider the operators for 

(Ofe, 3, a), and define the following operator on the bundle TM © (M x R): 

=X] /lfc2 (fe) J ( fe )^ u '^l £/ *)- 

k 

In view of Proposition 17. 2i P is a vfDO of order —1 with principal symbol 

Cndiag {w\-w>w\)- 

The operator A : C°°(TM) -> C°°(TM), 

A = -c„V(-A)- 3/2 <5, 
is a ^DO of order —1 with principal symbol 
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Therefore, the principal symbol of the operator 

P + diag(A, 0) 

equals 

c„dia g (%1-Mr 1 ), 

which means that P + diag(A, 0) is an elliptic ^DO of order —1. 
Now, the same arguments as in [29] [32] [31] show that the operator 

r M (P + diag(A,0)) : H S (U) -» H S+1 (M) 

has closed range and finite codimension. 
Since on H S (U) 

r M P = r M Zu I u, (7.5) 
(|7.3[) and the argument of GM OH] show that the equation 

r M (P + diag(A,0))[«,^] = [v,tp] 

has a solution [it, ^] € 7i s (U) for every pair [i/, <£>] £ H S+1 (M). Then by (|7.5|) 

[v^} + [Vf 1 0}=r M r u l u [u^] 

with / = c n (— A) _3 / 2 5u, which proves the lemma. 

Continuing the proof of Theorem l7.3[ observe that by Lemma l7.4( for every pair 
[v, if] S C°°{M) there exist [u, ip] G C°°(U) and / € C*°°(Af) such that 

[« + V/ > v>]=r J jf2^2i 7 [ti,# (7.6) 
For (x,0 € SZ7, define ^(x.f) by 

tt^ (a, 0) e flf, ^ < 0, ^+ > 0. 

Define 

™ + (z, £) = / («i(7x,t(*)7i, e (*) + ^(7x,«(<))) 
Jo 

«J~(z,0 = / («i(7»,e(*)'Vi,c(t)+V'(7x,c(t)))*- 

These functions belong to C co (SU). Let w = (w + + w~)\q + sm ■ It is easy to see 
that (w + + w~) is constant on the orbits of the magnetic flow; therefore, w* — 
(w + +w~)\sm and G C°°(SM). Since Tjj [u, ip] = (w + + w~)\q + su> we see from 
(|7.6|) that l*w = [v + V/, ip], which completes the proof of the theorem. □ 

7.2. Scattering relation and folds. The main aim of this section is to prove the 
characterization (|7.4[) for the space C^(d + SM) (Lemma l7.6p . We will proceed in 
the same way as in [32] . 

Preserving the notations of the previous section, define the map 

®:d(SM)^d-SU 

by 

9(x,0=^ M (.x,0, (x,0ed(SM). (7.7) 
Since £ + is smooth in SU, <I> is smooth as well. 

Lemma 7.5 (cf. [32] Theorem 4.1]). $ is a fold map with fold S(dM). 
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We recall that a smooth map / : M — > J\f between two smooth manifolds M. 
and M of the same dimension is said to be a Whitney fold (with fold L) at a point 
m G L if / drops rank by one simply at m, so that {x : df{x) is singular} is a 
smooth hypersurface near m and ker(<i/(m)) is transverse to T m L. 

Example. Let S be a manifold, 7Z a smooth function on S having as a regular 
value, and M. = |7?. _1 (0)}. Let X be a nonzero vector field on E such that 
XTZ(m) — 0, and XXlZ(m) / for a point m G A4. Let TV" be a hypersurface in £ 
transversal to X such that the projection / : A4 — > TV along the integral curves of 
X is well denned. Then / is a Whitney fold at m with fold M n {(XTe)" 1 ^)}. 

Indeed, we may assume without loss of generality that S is a domain in K n , 
X = d/dx n , and J\f is the hyperplane {x n — 0}. In this case, we easily verify the 
claim by straightforward calculations. 

Proof of Lemma \ 7. 5\ Let p be a defining function of M in U with gva,dp(x) = v(x) 
for x G 0M . Then 

G AI (/7O7r)0,£) = (£,gradp(:r)), 

G£(p 70(3,0 = (&V € gradp(x)} + (Y x (u), gradp(a;)}. 

Therefore, for (a;,£) G S(5M) we have G M (poir)(x,£) = 0, while G^(p o 7r)(x, £) = 
— A(x, + (Y^(^), v(x)) ^ by strict magnetic convexity of 9M. 

We arrive directly at the above Example if we take S = 5M , 7£ = p o 7T, A 7 ! = 
d(SM), M = d(SU), and A = G M . This completes the proof of the lemma. □ 

Define the extension operator 



where S is the scattering relation. 

Lemma 7.6 (cf. [321 Lemma 1.1]). If(M,g,a) is a simple magnetic system, then 



converse. If Aw G C°°(d(SM)), then from Lemma O and [13 Theorem C.4.4] 
we deduce the existence of a smooth function v on a neighborhood of the range 
$(d(SM)) such that w = v o $. 

Let <f : SM -> dSM denote the map 

*(x,0=^ ) (^0, (x,Z)€SM, 
and : 5t/ -> 9_SJ7 the map 

<Mz,0=^ +(x '« (*,£), M)gSC/. 

Note that iu" = o '!>. Therefore, to" =»o$o S^ 1 o vp. It is easy to see that 
$o5 _1 o* = ^u\ S M- Since is smooth on SM, we conclude that w s G C°°(SM), 
Le.,weC™(d + SM). □ 



A : C(d+SM) 



C(SM) 



by 




G 
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7.3. Hilbert transform. From now on, we let (M, g) be an oriented Riemannian 
surface. Given v € T X M, we will denote by v±_ the vector obtained by rotating v 
by 7r/2 according to the orientation of M. In coordinates (v±)i — £ijV J , where 

1 



Consider the fiberwise Hilbert transform [32, (1.4)] 

H < x ' ® = ^j s m 1 ^^ ) u(x, rj) da x ( V ), £ e S X M. 

If we fix x 6 M and a reference point a S S X M 1 any function on the fiber can 
be treated as a function of an angular variable. Then 

Hu(x, 9) = — / cot I — - — I u(a;, 99) d^. 



2^ 7o 
Define 

(G ± u)(a; > = (a,Vlti> = -<e > Vx«) J 
where V±u = e\/'u. The following commutation formula holds [321 Theorem 1.5]: 

[H,G]u = G ± (u ) + (G±u) , (7.8) 

where 

u o(x) = — u(x,£)da x (£) 

/7r JS.M 

is the average value on a fiber. 

Let V be the infinitesimal generator of the action of S 1 on the fibers of the 
canonical projection tt : SM — * A/. Then the generator of the magnetic flow -0* is 
given by 

G M = G + AF, 

where A is a function on M such that 17 = Af2 a , with f2 a the area form of g. 
Since H commutes with XV, we get from (I7.8|) 

[If, G> = Gx(w ) + (Gj.«) . (7.9) 

Substitute u = w e C™(d + SM), into (J72J). This yields 

G M W - -G ± (w»)o - (G ± w%. (7.10) 

Define the operator 

B : C(SM) C(d+SM) 

by 

Bu(x, = u(x, 0-uo S(x, £), (x, £) e d+SM. 

Clearly, 

JG^u = -Bu. (7.11) 
Using (|7.1ip . we deduce from (|7.10|) the identity 

BHAw = I(G ± (w%) + I{(G ± w%), (7.12) 

since w'|o(sm) = -Aw. 
Note that 
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and 

(G_lw b )o = -—S±Xlw, 

Z7T 

where 

S±v = — 

Hence, we can rewrite f|7. 1 2[) as 

BHAw = -^-l[V±I^w,S ± I^w}. (7.13) 

Z7T 

7.4. Dirichlet-to-Neumann map. Given a compact Riemannian manifold (M, g) 
with boundary, denote the Laplace-Beltrami operator associated with g by A g . 
Consider the Dirichlet problem 

A g u = on M, u| aM = /■ 

The DN map is defined by 

A g (f) = (v,Vu\ aM ). 

In the two dimensional case, the DN map can also be described as follows. Let 
(h, ft.*) be a pair of conjugate harmonic functions on M, 

Vh = Vj_h*, Vh* = -Vj_h. 

Let h° and h° denote their traces on DM. Then 

A(ft°) = (v,Vh% M ) = -{v,(Vxh)\ aM ) = K,(V/i)| aM ) = (v±,Va M h ), (7.14) 

where Vqm is the gradient w.r.t. the induced metric on dM . 

The following theorem is an analog of [32l Theorem 1.3] and states that the 
scattering relation of a simple magnetic surface (M, g, a) completely determines 
the DN map of the metric g. 

Theorem 7.7. Let (g,a) and (g',a r ) be simple magnetic systems on a compact 
surface M with boundary such that g\dM = g'\dM- Assume that the scattering 
relations S and S' of these systems coincide. Then A g = A g i . 

Proof. Assume that h, /i* is a pair of smooth conjugate harmonic functions on M. 
Then G±h = Gh* = G M /i*. By Theorem Ol there are w £ C^{d+SM) and 
/ S C°°(M) satisfying l^w = h, X\w = V/. From (j7TT3f we then obtain 

BHAw = ~Bh°„ (7.15) 

ZTT 

since 6±Vf = 0. Hence, the following holds: 

Lemma 7.8 (cf. [32|, Theorem 1.6]). If h,h r is a pair of smooth conjugate har- 
monic functions on M, then there is w 6 C^{d+SM) such that h — IqW and 
equation (|7.15jl holds with h® the trace of h* on dM. 

In the opposite direction we have: 

Lemma 7.9 (cf. [12 Theorem 1.6]). Suppose h° £ C°°(dM) and w £ C™(d+SM) 
satisfy equation (|7.15p . Define h := TqW and let be the harmonic continuation 
of hi to M. Then h and h* are conjugate harmonic functions. 
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Proof. Let q be an arbitrary smooth extension of h® to M. Note that 

Using (|7.13p , we can therefore rewrite (|7.15[) as 

-l[V ± h,<p]=l[Vq,0], 

with (p = 8±I\W- 
Thus, 

I[V±h + Vq,<p] = 0. 
By (17. 2|) we then have <p = and Vq + Vj_/i = Vp for some smooth function p on 
M with p\dM = 0. Therefore, /i and q — p are conjugate harmonic functions. Since 
(<Z — p)|sa/ = ft°,j we have q — p = /i*, which implies that /i and /i* are conjugate 
harmonic functions. □ 

Continuing the proof of Theorem [73 we have the following procedure to obtain 
the DN map from the scattering relation. For an arbitrary given smooth function 
h® on dM we find a function w G C^ D (d+SM) that solves equation (|7.15[) . Then 
the functions ft. = 2tt(Aw)o (notice that 2tt(Aw)o — I^w\qm) and h® are traces of 
conjugate harmonic functions. This gives the DN map by means of (|7.14[) . □ 

7.5. Proof of Theorem IT. X I We end up with the proof of this theorem in the 
same way as in |32j . If two simple magnetic systems (g, a) and (g' , a') on a compact 
surface with boundary have the same boundary action functions, then by Theorem 
I2.2l we may assume that g\dM = g'\dM and, by Lemma [2.51 their scattering relations 
coincide. Theorem 17.71 then tells us that the DN maps of the metrics g and g' 
coincide. Now, the result of [2lJ [19] implies the existence of a diffeomorphism 
/ : M — » M , which is the identity on dM, and of a function oj such that g 1 = w 2 f*g. 
Next, Theorem 16.11 yields lo — 1 and a' = f*a + dip for a smooth <p vanishing on 
dM. This concludes the proof of the theorem. 

Appendix A. Geometry of magnetic systems 

A.l. Mane's critical value and simplicity. Here we adapt a certain part of the 
theory of convex superlinear Lagrangians to the case of manifolds with boundary. 

Let M be a compact Riemannian manifold with boundary and let L : TM — > M 
be a C°° Lagrangian satisfying the following hypotheses: 

• Convexity: For all x £ M the restriction of L to T X M has everywhere 
positive definite Hessian. 

• Superlinear growth: 

L(x,v) 
lim — -- — = +oo 

\v\— too \v\ 

uniformly on x G M. 
The action of L on an absolutely continuous curve 7 : [a, 6] — > M is 

Az(7)= / L( 7 (t),j(t))dt. 

J a 

For each k G M., the Mane action potential A/- :MxM->lU {— 00} is defined 

by 

A k (x,y) = inf A L+fe (7), 
7eC(z,j/) 
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where C(x,y) = {7 : [0, T] —* M : 7(0) = x, 7(T) = y, 7 is absolutely continuous}. 
The critical level c — c(L) is defined as 

c(L) = sup{fc £ M : AL + k("f) < for some closed curve 7} 

= inf{£; £ E : AL+k(l) > for every closed curve 7}. 

Proposition A.l. For k < c(L), Ak(x,y) = —00 for all x,y £ M. For k > c(L), 
Ak(x, i/)el for all x, y £ M. 

Proof. The same as in [6j Proposition 2-1.1]. □ 

One more characterization of c(L) is useful. Recall that the Hamiltonian H : 
T*M — > R associated with L is defined by the Fenchel transform 

H(x,p) = sup{p(v) - L(x,v) : v £ T X M}, 

and the supremum is achieved at v such that p = ^(x,v). 

Proposition A. 2. // there exists a C 1 function f : M — * K such that H(df) < k, 
then k > c(L). 

Proof. The same as in Lemma 5]. □ 

Recall that the energy function E : TM — > M for L is defined by 

dL 

E(x, v) = — (x,v) ■ v - L(x, v), 
ov 

and that the energy function is constant on every solution x(t) of the Euler- 
Lagrange equation 

l^-( x (t),x(t)) = ^(x(t),x(t)). (A.l) 

Let ip* : TM -> TM be the Euler-Lagrange flow, defined by ip f (x, v) = (7^), j(t)), 
where 7 is the solution of (jA.ll) with 7(0) = x and 7(0) = v. For x £ M and 
k £ M, the exponential map at x of energy k is defined to be the partial map 
exp x : T X M — > M given by 

exp^(to) = 7to ip^v), t > 0, v £ T X M, E(x,v) = k. 

Then expj; is a C^^-smooth partial map on T X M which is C°°-smooth on T x M\{0}. 
The next proposition is similar to [8] Theorem D] and has a similar proof. 

Proposition A.3. // cxp£ : (exp£) _1 (M) -» M is a diffeomorphism for every 
x £ M, then k > c(L). 

Proof. Fix q £ M. Given x £ M, let ^i q _ x : [0, T q}X ] — » M be a solution of the 
Euler-Lagrange equation with energy k, joining q to x. Consider the function 
f{x) = Afc(7 9i2; ). It is easy to sec that the assumption of the lemma implies this 
function is smooth in M \ {q}. It follows from the first variation formula of [8, 
Lemma 4] that for x £ M \ {q} 

dL 

d x f(w) = -^(x,% tX (T qtX )) -w, 

which implies 

H(x,d x f) = E(x, A f qtX {T q:X )) = k. 

This last equation also implies that \d x f\ is uniformly bounded for all x £ M \ {q}. 
Thus / is Lipschitz in M and a smoothing argument as in Q3] shows that for 
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any e > there exists / £ C ca (M) for which H(x,d x f) < k + e for all x £ M . 
Thus by Proposition IA2] k > c(L). □ 

The next proposition is an analog of [H Proposition 3-5.1] and has the same 
proof. 

Proposition A. 4. If k > c(L) and x,y £ M x ^ y, then there is 7 £ C(x,y) such 
that 

A k (x,y) = A L+fe ( 7 ). 
Moreover, the energy of 7 is £(7,7) = k. 

Now, we apply the above to the case of magnetic systems. For a simple magnetic 
system (M, g,a), the magnetic flow can also obtained as the Euler-Lagrange flow 
with the corresponding Lagrangian defined by 

L(x,v) = -\v\ 2 g - a x (v). 

Lemma A. 5. Let (g, a) be a simple magnetic system on M . For x, y £ M , x ^ y, 



&1/2(X,V) = ^L+l/2{lx,y) = T Xty - 




where j Xt y : [0, T x ,y] —> M is the unit speed magnetic geodesic from x to y. 

Proof. It is easy to see that the simplicity assumption implies that for this La- 
grangian the assumptions of Proposition [A~3l hold for all k sufficiently close to 1/2. 
Therefore, the proposition gives 1/2 > c(L). Then Proposition IA.4I shows that, 
given x y in M, there is 7 £ C(x,y) with energy 1/2 (i.e., 7 is parametrized by 
arc length) such that A(x,y) = A( 7 ). Using simplicity, one can then prove that 7 
is a unit speed magnetic geodesic, i.e., 7 = 7 x ,j/- Q 

A. 2. Magnetic convexity. Let M be a compact manifold with boundary, en- 
dowed with a Ricmannian metric g and a closed 2-form fi. Consider a manifold 
Mi such that M[ nt D M. Extend g and fl to Mi smoothly, preserving the former 
notation for extensions. We say that M is magnetic convex at x £ dM if there 
is a neighborhood U of x in Mi such that all unit speed magnetic geodesies in U, 
passing through x and tangent to dM at x, lie in Mi \ M mt . If, in addition, these 
geodesies do not intersect M except for x, we say that M is strictly magnetic convex 
at x. It is not hard to show that these definitions depend neither on the choice of 
Mi nor on the way we extend g and to Mi. 

As before, we let A denote the second fundamental form of dM and v{x) the 
inward unit normal to dM at x. 

Lemma A. 6. If M is magnetic convex at x £ dM. then 

A(x, v) > (Y x (v),v(x)) for all v £ S x {dM). (A.2) 

// the inequality is strict, then M is strictly magnetic convex at x. 

Proof. Suppose M is convex at x. Choosing a smaller U if necessary, we may assume 
that there is a smooth function p<mU such that | gradpj = 1 and dMflU = p _1 (0). 
Further we may assume that all the above geodesies lie in U~ = {x : p(x) < 0}. 
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Let v G S x (dM) and 7(i) be the magnetic geodesic with 7(0) = x, 7(0) = v. By 
assumption, p o j(t) < for all small t. Therefore, 

d 2 



< 0. 

t=o 



Since 

^[P°7(*)] =^(gradp(7(«)),7(*)> 

= (V i(t) gradp(7(*)),7(t)> + (gradp(7(«)),7(*)> 

= Hess 7(t) p(7(t),7(t)) + (gradp( 7 (t)), F( 7 (t))) 

and since A(ai, £) = — Hess^ p(v, v) and gradp(.x) = v(x) when £ d(SM), we 

arrive at (|A.2|) . 

Now, assume that (|A.2[) is strict. Then there is 6 > such that for every 
magnetic geodesic 7 in TV with 7(0) = x and 7(0) = v G S x (dM), 



t=o 



Thus, there is a small e > such that 

po 1 (t) <~^t 2 for alU G (-£,£). 
This implies the second statement. □ 
A. 3. Exponential map. 

Lemma A.7. Trae map exp£ : T X M -> M is C 1 and C°° on T x M\{0}. This map 
is C 2 if and only if f2 = 0. 

Proof. Recall that exp^(iw) = it o ip f (x,v) :— j(t,x,v) where t; has norm one and 
t > 0. Clearly this implies that exp£ is C°° on T X M \ {0}. Since 



d 
~dl 



7T O 1p (X, V ) = V 



we see that the directional derivative of exp^ at in the direction of v exists and 
equals u. In a coordinate system around x write 

exp^(to) =Y{t,x,v) (A.3) 

where j l depends smoothly in (t, x, v). Differentiating with respect to v we obtain 

^^(tv) t =p:(t,x,v). 

Since 7 l (0, x, v) = x l , §^j(0, x, v) = 0. Therefore, 

c^exp^' 4 1 0^f l 

lim — - x (tv) — lim - -— (t, x, v) 

1 ^,x,v) = -J-(0,x,v)=^ 



dtdvi ' ' dv^dt 

since 7*(0, x, v) — v 1 . Since S X M is compact, the above limit is uniform in v G S X M 
and thus exp^ has continuous partial derivatives at 0, i.e., exp£ is C 1 . 
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Suppose now that exp£ is C 2 . Differentiate (|A.3|) twice with respect to t (v £ 
S X M) 

ov J v 

Using the equations of a magnetic geodesic 

■■i i • 7 • k-v^i • k 

7 +T1 T jk = Y kl 

we obtain 
Let t -> 0+. Then 

..i„.* /" d 2 exp£' 



-(0) =y fc i ( ! r)t;*. 



Since this holds for all v £ S X M we must have Y£(x) = 0, i.e., f2 = 0. □ 

A. 4. Santalo's formula. If (M, g) is a compact Riemannian manifold with bound- 
ary, we endow its unit sphere bundle SM with its usual Liouville (local product) 
measure dS 2rt_1 , and endow the bundle d+SM with its standard measure dT, 2n ~ 2 
(again a local product measure where the measure on the fiber is the measure of a 
hemisphere). Denote by dfi the measure on d+SM given by 

d»(x,0 = (e,^))ds 2 "- 2 (x,o, 

where v(x) is the inward unit normal of dM at a point x. 

The following version of Santalo's formula holds for magnetic flows. 

Lemma A. 8. Suppose that (M, g, a) is simple. Then for every continuous function 
ip : SM — > R we have 

[ <pdY? n - l =( dfx(x,0 f"' ip{ lx<i {t),% A {t))dt. (A.4) 

JSM Jd+SM JO 

Proof. The argument we use is the same as in [34]. We give it for the sake of 
completeness. 

First, we recall the well-known fact that the Liouville measure is invariant under 
the magnetic flow (for example, because ajQ An = (ujq + 7r*f2) A ™ while ujq + 7r*£l is 
flow invariant). Now, let D = {(x,£;t) £ d+SM : < t < £(x,£)}, and define 
* : D — > SM by ^(x, £; t) = ip l {x, £), where ip l is the magnetic flow. Then 

f ipdY? 11 - 1 = [ ((po*)**(dE 2n - 1 ). (A.5) 

JSM JD 

By construction, if? conjugates ip 1 with the flow generated by d/dt on D. Since 
(il] 2 " -1 is invariant under ip*, ^*(<iE 2 " _1 ) is invariant under the flow of d/dt. Then 

**(dS 2n_1 )(a;,^t) = a(x,C)dS 2 "^ 2 Adt 

for some function a on d+SM, so that (|A.5[) takes the form 

r /•<(*.€) 
<pdX 2n - x = a(x,0d^ 2n ' 2 (x,0 P(7*,c(*).-V*,«(*))*- 

SM Jd+SM Jo 

We are left with proving that a{x, £) = (£, v(x)). To this end, it suffices to show 
that 

^{dJ? n - x ){x, £; 0) = (e, i/(a:)> d£ 2n ~ 2 A dfc. 
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Note that dS 2 ™- 1 = -dS 2 ™" 2 A dr on d(SM), where r(x) = dist(x,<9M) is 
the distance function from x to dM. Since the differential of ^ is the identity on 
T( x .£)(d+SM) and takes d/dt to the generator G^ of the flow ip*, we have 

**(dE 2n - 1 )(a;,^0) = (G^r) g?£ 2 "~ 2 (x, £) A dt. 
As soon as G^r = (£, Vr) = — (£, f(x)), we are done. □ 

A. 5. Index form of a magnetic geodesic. Let (M,g,a) be a simple magnetic 
system. For every x G M, exp£ : T X M — > M is a diffeomorphism restricted to a 
suitable set in T X M which is diffeomorphic to a closed ball. 

Let 7r : SM — ► M be the canonical projection and let for v G SM, 

V(v) :— ker d„7T, 

which is an (n — l)-dimensional subspace of T V SM, and 

E(v) := V(v) © iG (l («). (A.6) 

Lemma A. 9. If 7 : [0, T] —>■ M is a unit speed magnetic geodesic, then 

d m ^(E)nv(j(t)) = {0} 

for every t G (0, T]. 

Proof. Take w G SM and i G (0, T]. From the definition of exp^ one sees right 
away that 

image(d to exp£) = d^^-n:^^ 1 (E)). 
Since d w exp£ is a linear isomorphism for every w G T X M at which exp^ is defined, 
the lemma follows. □ 

Given a unit speed magnetic geodesic 7 : [0, T] — > M, let .4 and C be the 
operators on smooth vector fields along 7 defined by 

A(Z) = Z + R(j, ZYf - Y(Z) - (V Z F)( 7 ), 

c(z) = i?(7,z) 7 -r(z)-(Vzr)( 7 ). 

A vector field J along 7 is said to be a magnetic Jacobi field if it satisfies the 
equation 

A(J) = 0. (A.7) 
Let A denote the M-vector space of smooth vector fields Z along 7 such that 
Z(0) = Z(T) = 0. Define the quadratic form Ind : A -> M by 

Ind(Z, Z) = £ { |Z| 2 - (C(Z), Z) - (F( 7 ), Z> 2 } dt. 

Note that 

Ind(Z, Z) = - / {(.4(Z), Z) + <y(7), Z) 2 } di. 
Jo 

Clearly, Ind(Z, Z) generalizes the index form of a geodesic in a Ricmannian 
manifold. The next lemma is an analog of a well-known assertion from Ricmannian 
geometry. 

Lemma A. 10 (Index Lemma). If Z G A is orthogonal to 7, iften 

Ind(Z, Z) > 0, 
wif/i equality if and only if Z vanishes. 
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Proof. Note that the subbundlc E defined by (|A.6j) is Lagrangian. 

If £ £ E(v), then Jf (t) = dir o dip*^) satisfies the Jacobi equation (|A.7j) . Since 

d j(t)-*\ Em)) ■ E(i(t)) -» T 7(t) M 

is an isomorphism for all t £ (0, T], there exists a basis {£i, . . . , £„} for -E(i>) such 
that {%(*),..., Jg n (t)} is a basis of T 7 ( t) M for all t e (0,T]. Without loss of 
generality we may assume that £i = G M (n), = 7 and ./^(O) = for i > 2. 
Let us set for brevity Ji = . Then if Z is an element of A, we can write for 

te(o,T] 

n 
i=l 

for some smooth functions /1, . . . , /„. The functions fi can in fact be smoothly 
extended to t = 0. Indeed, for i > 2, we can write Ji(t) — tAi(t) where Ai is a 
smooth vector field such that Ai(0) = Ji(0). Since {j(t), A 2 (t), . . . , A n (t)} is now 
a basis for all t £ [0, T], there exist smooth functions <?i such that for all t G [0, T] 

n 

Z(f)=Si(*)7(*)+X)Si(t)^(t)- 

i=2 

Therefore for t E (0,T], gi(t) = h(t) and for i > 2, 9l (t) =tfi(t). Since Z(0) = 0, 
gi(Q) — for all i and the fi's smoothly extend to t = 0. 
Now we can write 

lnd(Z,Z) = -J2 f (AifiJ^fjJj) dt- [ (Y(j),Z) 2 dt. (A.8) 

An easy computation shows that 

A(fiJi) = fiJi + 2/4 - fiY(Ji) + fiA{Ji). 
Since Ji satisfies equation (|A.7[) , then A(Ji) = and hence, 

(AifiJi), J,) = (Ji, Jj) + 2fi{Ji, Jj) - fi(Y{Ji), Jj). 
Observe that since E is a Lagrangian subspace, 

{Ji: Jj) — (Ji, Jj) + (Y(Ji), Jj) = 0, 

and then 



(A(f i J i ),J j ) = j t (fi(Ji,J j )). 



Now we can write 



T 



(AtfiJjJjJj) dt= (fiJufjJj) - (f i J i ,f j Jj)dt 







T 







T 



dt. 



Combining the last equality with (|A.8[) we obtain 

Ind(Z,Z)= [ IVM 2 dt-(V/,4^) T - / (Y{i),Z? 
But Z(Q) = Z(T) = 0, therefore 

i(v,v)= f T \\y f,J, di- I in: >.zr,ii. iA-n 



Jo 11 i=1 Jo 
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Now let 

n 
i=2 

Since J\ = 7 we have: 

n n 

( E fi J *< E /* J *) = (/n + ^ A7 + w) = il + 2/1 <7, w> + (w, w). 

i=l i=l 

Differentiating (Z, 7} = we get 

(Z,j) + (Z,Y(j)) = 0. 

But 

n 

(Z,7> = (5^/^,7) = A + (W, 7) 

i=l 

since (Ji,7) = for all i. Therefore 

(r( 7 ),z) 2 = A 2 + 2/1 (w, 7) + (^, 7 } 2 - 

Thus 

n n 

( E /« J *. E - < y (7). z ? = (W, W) - (W, i) 2 . 

i=l i=l 

If we let II^ 1 - be the orthogonal projection of Wtcy 1 , the last equation and (|A.9|) 
give: 

lnd(Z,Z) = f \\W x fdt > 
Jo 

with equality if and only if vanishes identically. But if vanishes, then 

n 

-(^,7>7 + E/^ =0 

»=2 

which implies that the functions are constant for i > 2. Thus Z is of the form 
/i7 + •/ where J is a magnetic Jacobi field. But Z{T) = implies J(T) = 0. Since 
the J[s are linearly independent at T, J must vanish identically and since Z is 
orthogonal to 7, Z must also vanish. □ 

Appendix B. Study of a certain class of integral operators with 

SINGULAR KERNELS 

As we mentioned before, the fact that the magnetic exponential map is smooth 
in polar coordinates only forces us to work in polar coordinates as well. In this 
appendix we study a class of operators that naturally arise in our analysis. 

Let U C M™ be open and g be a smooth Riemannian metric in a neighborhood 
of U. 

Lemma B.l. Let A : Cq(U) — > C(U) be the operator 

Af{x)= I [ A{x,r,u;)f(x + rLu)drda x {u;), (B.l) 
Js x u Js. 
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with A e C°°{U x E x S X U). Then A is a classical ^DO of order -1 with full 
symbol 

oo 

o 

w/iere 



s x t/ 



Proof. Notice first that if A is an odd function of (r, u>), then .A/ = 0. Therefore, 
we can replace A above by A cvcn (r,ui) — (A(r,uj) + A(— r, —u>))/2. Next, it is easy 
to check that we can integrate over r > only and double the result. Therefore, 



Af{x) = 2 / A cvm {x,r,uj)f{x + ruj)drda x {uj), (B.2) 

JSxM JO 

Consider now r, oj as polar coordinates for z = ru>, and make also the change of 
variables y = x + z to get 

Af{x) = 2{detg{x)) 1 ' 2 f Acvcn (x, \y - x\ g , , d V> ( R3 ) 

J V \y- x \ g J \y - Ag 

where the subscript g refers to g(x). Let 

N-l 

A cvcn (x, r,uj) = 2J A eveaik (x, Lo)r k + r N R N (x, r, us) (B.4) 

k=0 

be a finite Taylor expansion of j4 cvcn in r near r = with N > 0. It follows easily 
that 2A even ^ k (x, (J) = Afe(a;,w) + (— l) h A k (x, — ui), where k\A k = d k \ r= oA, and in 
particular, A cvcn . k (x, cu)r k is even w.r.t. (r, (S). The remainder term contributes to 
(|B.3|) an operator that maps L^. omp (U) into H N ~ N °(U) with some fixed ATq. To 
study the contribution of the other terms, write 

•4even,fc/(z) - 2(det g(x))^ 2 J A even , fc (x, ) [y - x\ k g - n+1 f(y) dy. (B.5) 

The kernel of -4 C von,fc is therefore a function of x and z — y — x, with a polynomial 
singularity at y — x = 0, and it is therefore a formal ^PDO with symbol that can be 
obtained by taking Fourier transform in the z variable. Motivated by this, apply 
the Plancherel theorem to the integral above to get 

A cvcn , k f(x) = (2tt)-" J e***a k (x,Qfe)dS, 

where 

a k (x,0 = 2 J e -^<A cvc ^ k ^x,^^\y - x^+^detgix)) 1 / 2 dy 



S^U JO 



oc 



g - 

3 - l ™-«A cvcnifc (x,w)r fc drda x {u) 
' lrul< A k {x,uj)r k drda x (uj) 



S X U J —oo 



2m k / A k {x,uj)5 {k) {uj-C)da x {uj). (B.6) 

JS X U 
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In the third line, we used the fact that A evcllt k(x, uj)r k is even. Note that ak(x, £) is 
homogeneous in £ of order — k — 1 and smooth away from £ = but a distribution 
(in S') near zero. To deal with this, choose \ G Co° supported in |£| < 1 and equal 
to 1 near £ = 0. Write a(x, £) = x(O a ( iE >C) + (1 — XiO) a ( x '0- The second term 
is a classical amplitude, while the first one contributes the term 

A cvcn ,k(x*f) (B.7) 

to (|B.5|) that is smooth, as can be easily seen by making the change of variables 
z = y — x in (|B.5[) , see also (|B.16[) . □ 

Remark B.2. If A(x,r,uj) and g are smooth of class C k only, then A is an WDO 
with an amplitude of finite smoothness l(k), admitting a finite expansion. If k 3> 1, 
then 1^1, and one can still construct a finite order parametrix of an elliptic WDO 
in this class and the usual H Sl — > if 82 estimates still hold, if fc 3> 1, depending on 
si, S2- This has been used already in [571155] . 

We return to the analysis of the singular operator A introduced in Lemma IB. II 
under the assumption that A and g are analytic. Our reference for the calculus of 
analytic <I>DOs is [ID] . 

Lemma B.3. Let A : Co(U) — > C(U) be the operator W. 1\) with A(x,r,oj) analytic 
for (x,u>) G U x S X U , and r £ M smc/i f/iai cc + raj G {/. T/ien A is an analytic 
DO of order —1 with a symbol expansion as in Lemma \B.l[ 

Proof. Notice first that A is analytic pseudolocal, see [40] Theorem V.2.1]. 

By performing the change of variables ui' = g^/ 2 (x)LJ, we reduce the lemma to 
the case where g is the Euclidean metric. Let U' CC U. 

Let us estimate a/c(a;,£), see Lemma IB. II Since A(x,r,oj) is analytic, and is 
homogeneous of order — k — 1, we have 

\a k (x,0\ <C fc+1 fc!|£|- fc - x 

with some C > 0. Using the homogeneity, we get 

\dga k (x,£)\ < C fc+ l a l +1 a!fc!|£|- fc -l a l- 1 (B.8) 

for £ is in a complex neighborhood of K and x in a complex neighborhood of U' . 
Therefore, there exists a pseudoanalytic symbol a ~ ^ at, see [40] . This symbol is 
defined by 

oo 

ofoO = £v*(Oafc(z.O. ( B -9) 

fe=0 

where ifk have the properties (see [30] V.2]): < 9?^ < 1, <Ph(Q — f° r I CI < 
2i?max(fc,l), (p fc (£) = 1 for |£| > 3i?max(fc,l), |£> Q (^ fc | < (C7/J?)I Q I for |a| < 2k, 
where R > 1 is a large parameter. We will prove next that a(x, D) differs from A 
by an analytic regularizing operator. 

Let u G £'(U'). Let s > be such that u can be represented as a finite sum of 
derivatives of continuous functions of order not exceeding s. 

By pTIj) . for any AT > 1, 

N-l 

Au= ^2 a k (x,D)u + TZ N u, (B.10) 

fc=0 



GO 
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where, despite the strong singularity of au at £ = 0, a,k{x,D) have regular (inte- 
grable) Schwartz kernels, and 

TZnu(x) — 2 / / r N RN(x,r,Lu)u(x + ruj) dr da{uj) 
JS"- 1 Jo 

= 2 J \x-y\ N - n+x R N (x,\x-y\,^^)u(y)dy. (B.ll) 
We express Rn in its Cauchy form as 



1 

H N {x,r,uj) = jj-j: 

We have 



R N {x,r,uj) = —A J d?A cvcn ,k(x,tr, w)(l - t) N 1 dt. 



\D°"R, N u\ < C N a\ in U' for |a| < N - s. (B.12) 
Splitting the sum (|B.9p into two parts, we write 

JV-l oo 

a(x,D) = Op( J2 MOak(x,0)+Op( E MOa k (x,0)- (B.13) 

fe=0 fc=JV 

For the second term we have (see (3.15) in chapter V in (40] ) 

oo 

D a Op(Y2<p k (£)a k (x,Q)u<C N a\ in U' for |a| < N - s. (B.14) 

k=N 

We are left to compare the first sum in (|B.13j) with the sum in (|B.10j) : 

JV-l 

B N u := Op( i 1 - fk(0)ak(x,0)u= B' N u + B'^u, 
k=a 

where 

JV-l 

b' n u = o P ( (MO - ¥>*(0)«*(*> o' 

k=0 
JV-l 

S> - Op( 5^(1 - M0K(*>0)«- 

fe=0 

On supp (tfN — <fk)i we have 2Rk < |£| < 3RN provided that k < N and, as always, 
we assume that R 3> 1 . Using this and ()B.8|) , we get 

\D a B' N u\ < C(CRN)M- 1+S in U' for \a\ < N - 1, (B.15) 

compare with (3.17) in chapter V in [417] , 
We write B'^u in the form (see (|B.7|) ) 
JV-l . 

=2^2 4vcni (x, pr) (1 - <p N )'* f(z + x) dz. (B.16) 

This implies 

|£> Q £/>| < C^ORiV) 1 " 1 ^ in U' for |a| < JV - 1. (B.17) 
Combining (|Bl2|) , (|BTT4|) . [fRlS]) . and (|Bl7|) . we get 

|L> a („4 - a(x, D))u| < C N Nl in [/' for |a| < N - s. 
For N > s, choose \a\ — N — s to conclude that (A — a(x, D))u is analytic in U'. □ 
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